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ABSTRACT 
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This  document  contains  information  on  the  research  accomplished  under 
AFOSR  Grant  Ho.  A70SR  82^0089  during  the  time  period  1  February  1982 
through  31  January  1  983*  The  work  covers  several  different  areas  of 
optical  computing,  as  well  as  some  work  on  digital  processing  of 
optically  obtained  images.  The  primary  emphasis  of  the  work  is  on  the 
possible  applications  of  optics  to  interconnections  in  integrated 
circuit  technology.  Other  areas  of  effort  include  the  diagonalisation 
and  inversion  of  circulant  matrices  using  coherent  optics,  the  division 
of  complex  wavefronts  using  four-wave  mixing,  and  the  suppression  of 
speckle  in  coherently  formed  images.  Publications  during  the  last  year 
arising  out  of  work  supported  by  the  grant  are  also  detailed. 
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I.  INTRODUCTION 


This  report  covers  the  work  performed  on  AFOSR  grant  No.  82-0089 
during  the  time  period  1  February  1981  through  31  January  1982.  It  is 
divided  into  six  sections,  the  first  of  which  is  this  Introduction. 
Immediately  following  we  summarise  our  current  status  of  the  project 
aimed  at  investigating  the  applications  of  optics  to  the  interconnection 
problem  in  integrated  circuit  technology.  Section  III  deals  with  the 
results  of  our  work  on  the  use  of  coherent  optical  processing  for  the 
diagonalization  and  inversion  of  circulant  matrices.  Section  IV 
discusses  a  new  project  started  this  year,  namely  one  concerned  with  the 
possibility  of  performing  complex  wavefront  division  using  four-wave 
mixing.  Section  V  reports  the  early  results  of  another  new  project,  one 
aimed  at  suppressing  speckle  in  coherently  formed  images  by  means  of 
digital  image  processing  techniques.  Finally,  section  VI  details  the 
publications  and  meeting  presentations  arising  out  of  work  supported  by 
the  grant  during  the  past  grant  year. 
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II.  OPTICAL  IFTERCOHNECTIOHS  FOR  VLSI 


As  integrated  circuit  technology  advances,  device  sizes  are  being 
reduced  at  the  sane  tine  that  overall  chip  sizes  increase.  The  increased 
complexity  of  such  chips  places  ever  increasing  denands  on  interconnect 
technology.  The  RC  tine  delays  and  IB  voltage  drops  of  conventional 
aluninun  and  polycrystalline  silicon  interconnections  nay  ultinately 
eclipse  the  gains  anticipated  in  device  perfornance. 

This  project  is  concerned  vith  the  possibility  of  using  optical 
inaging  interconnections  to  overcone  sone  of  the  problens  inposed  by 
current  interconnect  technology.  The  substitution  of  optical  waveguides 
for  netallic  interconnection  lines  could  potentially  reduce  the  effects 
of  BC  tine  delays  and  IB  voltage  drops  now  experienced.  However,  our 
goal  is  to  evaluate  the  feasibility  of  using  optical  interconnects  that 
rise  above  the  chip,  rather  than  sinply  using  waveguides  to  ninick  the 
quasi-planar  topology  of  conventional  interconnections.  Ve  envision  the 
use  of  holographic  optical  elenents,  either  transnissive  or  reflective, 
to  provide  a  nultitude  of  non- interfering,  efficient  parallel  connunica- 
tion  channels,  either  into  a  chip  from  the  outside  world  or  between  two 
chips.  Eventually  the  problem  of  providing  interconnects  within  a  single 
chip  would  also  be  of  interest.  Since  the  particular  interconnect 
pattern  achieved  is  determined  by  the  holographic  optical  element  used, 
there  exists  the  tantalizing  possibility  of  dynamically  changing  the 
interconnect  pattern  to  meet  the  current  needs  of  the  processor 
performing  the  computations  on  the  chip. 

Of  the  various  problems  that  could  be  tackled,  that  of  communicating 
by  means  of  optics  from  the  outside  world  into  a  chip  appears  to  be  the 
most  amenable  to  current  technological  solution.  For  such  problems, 
there  need  be  no  flow  of  optical  information  out  of  the  chip,  and  hence 
the  difficult  problem  of  providing  on-chip  optical  sources  or  modulators 
is  avoided.  Considerable  time  was  wasted  in  the  early  part  of  this 
project  in  attempting  to  find  ways  to  integrate  small  liquid  crystal 
light  modulators  on  a  silicon  chip.  It  was  finally  concluded  that, 
while  the  goal  was  probably  achievable  with  considerable  work,  the  speed 


of  the  resulting  modulators  vould  be  so  slow  that  there  would  be  little 
interest  in  the  end  result. 

Providing  detectors  on  the  chip  poses  no  fundamental  problem.  As  the 
complexity  of  chips  grows,  the  needs  for  connections  into  the  chip  also 
grows.  Yet  current  pin- connecting  approaches  are  generally  limited  to 
providing  on  the  order  of  100  connections  to  the  outside  world.  The 
success  of  optics  in  meeting  these  needs  will  depend  critically  on  the 
sice  and  performance  of  the  detectors  that  can  be  realised  on  the  chip, 
and  on  the  stringent  positioning  tolerances  between  the  imaging  element 
and  the  chip. 

Our  work  on  the  imaging  interconnect  problem  has  three  separate 
aspects.  One  deals  with  the  optical  issues  associated  with  the 
interconnect  devices  (holograms),  and  the  second  with  the  on-chip  photo¬ 
detectors  that  provide  the  means  for  inputting  data.  The  third  consists 
of  a  continuing  search  for  algorithms  and  specific  problems  that  demand 
input  of  large  amounts  of  data  in  parallel  (while  requiring  rather 
smaller  output  capability).  A  fourth  aspect,  that  of  finding  ways  to 
place  optical  sources  on  chips,  is  viewed  as  a  longer-range  goal  and  is 
under  study  in  the  Integrated  Circuits  Laboratory  under  separate 
support.  It  should  be  mentioned  that  some  experts  in  this  field  feel 
that  in  the  long  term  a  solution  to  this  difficult  device  problem  will 
be  found,  possibly  by  using  a  buffer  material  (such  as  germanium) 
between  silicon  and  gallium  arsenide  to  provide  adequate  matching  of  the 
lattice  dimensions. 

It  is  attractive,  for  reasons  of  cost  and  simplicity,  to  consider 
light-emitting  diodes  (LED's)  as  prime  candidates  for  the  modulated 
sources  that  will  be  used  for  the  optical  communication  channels  into 
the  chip.  However,  the  spectral  purity  of  LED's  is  far  poorer  than  that 
of  lasing  diodes,  and  therefore  some  understanding  of  the  resolution  of 
holographic  optical  elements  when  used  with  LED  sources  must  be  gained. 
As  a  goal  towards  this  end,  a  system  has  been  under  construction  which 
we  hope  will  provide  a  test  bed  for  answering  questions  of  this  kind. 
Also  to  be  gained  from  these  experiments  is  knowledge  about  the 
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scattering  of  light  by  the  holographic  optical  elements  and  by  the 
silicon  chip  itself.  A  number  of  GaAsP  unbonded  LED's  were  obtained  from 
Hewlett  Packard,  and  mounted  in  an  array  format  on  a  DIP  package. 
Several  encapsulating  materials  are  being  examined  to  determine  which 
will  minimise  light  scattering  and  beam  distortion  and  provide  environ¬ 
mental  protection.  Hewlett  Packard  also  supplied  a  number  of  unbonded 
silicon  photodetec tors  of  the  type  used  in  electro-optic  couplers.  These 
have  also  been  mounted  in  DIP  packages  and  will  be  calibrated  and  used 
to  test  S/H  and  response  threshold  levels.  A  Fairchild  I-scan  CCD  line 
scanner  will  be  used  to  monitor  images  in  the  hologram  image  plane. 
Initially,  transmission  and  reflection  type  holograms  will  be  prepared 
using  conventional  sources,  geometries  and  films.  Later,  experiments 
will  be  conducted  with  image  plane  holographic  geometries  and  dichro- 
mated  gelatin  recording  materials.  These  systems  have  promising  white- 
light  reconstruction  and  high  diffraction  efficiency  characteristics. 

Preparations  for  these  experiments  have  now  been  largely  completed 
and  the  experiments  themselves  are  anticipated  to  be  under  way  in  the 
very  near  future.  In  the  meantime  the  problem  of  integrating 
photodetectors  onto  the  chip  will  be  under  study  in  parallel.  Initial 
efforts  in  this  regard  will  be  aimed  at  studying  fundamental  limits  to 
the  sise,  speed  and  sensitivity  of  these  devices.  Our  approach  is  to 
identify  the  problems  that  are  anticipated  when  attempting  communication 
into  a  chip,  and  to  gather  the  information  necessary  to  assess  the 
seriousness  of  these  problems  in  a  methodical  way.  We  intend  to  focus 
attention  on  fundamental  issues,  such  as  the  radiometry  of  the*  problem, 
which  determines  the  available  light  power  at  the  detectors,  and  the 
fundamental  limits  to  detector  sise  and  sensitivity. 

Finally  we  mention  a  new  concept  that  has  arisen  during  the  past 
grant  year.  Consider  a  single  source  of  information  that  must  be 
supplied  to  a  number  of  sites  on  a  chip.  Further  suppose  that  the 
algorithm  being  realised  on  the  chip  requires  that  bursts  of  information 
be  sent  to  different  sites  on  the  chip  (i.e.  with  different  interconnect 
patterns)  as  a  function  of  time.  Such  requirements  can  be  met  in  the 
following  way.  A  single  high-speed  optical  channel  is  to  be  fed  into 
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the  chip,  but  it  should  be  interconnected  to  different  sets  of  locations 
on  the  chip  as  a  function  of  tine.  We  can  satisfy  such  a  requirement  by 
use  of  an  acousto-optic  bean  deflector,  and  a  two-dimensional  array  of 
holographic  optical  elements,  each  providing  a  different  interconnect 
pattern.  The  hologram  functions  much  as  a  holographic  memory.  To 
address  a  particular  set  of  desired  points  on  the  chip,  the  acousto- 
o-ptic  deflector  sends  the  beam  to  the  particular  location  on  the 
hologram;  the  image  read  out  is  a  series  of  spots  at  the  desired  loca¬ 
tions  bn  the  chip,  where  detectors  convert  the  burst  of  high-speed 
optical  modulation  into  electrical  signals.  The  beam  deflector  then 
sends  the  light  to  another  hologram  in  the  array,  which  then  sends  a  new 
set  of  light  spots  to  presumably  another  set  of  detectors  on  the  chip, 
where  the  new  burst  of  information  is  converted  to  electrical  signals. 

The  concept  described  above  appears  to  be  a  powerful  one,  but  the 
missing  link  at  present  is  the  lack  of  any  well-defined  problem  or 
algorithm  that  requires  this  type  of  communication.  Nonetheless,  we 
feel  confident  that  such  algorithms  can  be  found,  and  we  will  be 
searching  for  examples  during  the  coming  grant  year. 


III.  COHEBENT  OPTICAL  TECHNIQUES  FOB  DIACONALIZATION  AND  INVEBSION  OP 

CIRCULANT  MATRICES. 

For  the  past  year  we  have  been  studying  the  possible  use  of  coherent 
optical  systems  for  diagonal! zation  and  inversion  of  circulant  matrices. 
Circulant  matrices  are  those  for  which  each  successive  row  is  a  simple 
circular  shift  of  the  row  above  by  a  single  element.  For  example,  in 
the  matrix  A  below,  the  numbers  1,2,3,  and  4  stand  for  four  distinct 
elements;  the  organization  of  those  elements  in  the  circulant  matrix  is 
as  follows: 
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12  3  4 
4  12  3 
3  4  12 


(D 


2  3  4  1 

A  remarkable  property  of  circulant  matrices  is  that  they  are 
diagonalized  by  the  discrete  Fourier  transform  (DPT),  the  resulting 
diagonal  elements  being  the  complex  eigenvalues  of  the  original  matrix. 
Thus  if  the  DFT  matrix  V  is  defined  by  (again  illustrating  with  a  4x4 
example) 

1111 
1  w  w2  w^ 

(2) 

1  w2  w^  w® 

1  w^  w^  w^ 

where  w  *  exp(-i2  tt/N),  we  have 

0  0  0 

A  «  V1  A  V  «  0  X2  0  0  (3) 

0  0  x3  0 
0  0  0  x4 


t 


« 


where  the  X  's  are  the  eigenvalues  of  A. 

Before  embarking  on  a  more  detailed  description  of  the  work,  it  is 
perhaps  worthwhile  addressing  the  question  why  would  anyone  want  to 
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diagonalize  circulant  matrices,  and  why  might  coherent  optics  be  a 
useful  approach  in  some  cases.  A  great  many  signal  processing  problems 
require  the  inversion  of  Toplitz  correlation  matrices.  Toplitz  matrices 
have  less  structure  than  circulant  matrices;  they  simply  have  the 
property  that  each  subdiagonal  of  the  matrix  has  elements  that  are 
identical.  Such  matrices  must  be  inverted  in  order  to  determine  the 
structure  of  optimal  linear  filters,  and  such  inversions  must  be  per¬ 
formed  frequently  in  environments  where  the  statistics  are  changing  with 
time  or  are  gradually  becoming  known  as  time  progresses.  Usually  the 
Toplitz  matrices  of  concern  are  banded  matrices;  that  is  only  the  main 
diagonal  and  a  few  subdiagonals  have  significant  non-zero  value.  It  can 
be  shown  that  a  large  banded  Toplitz  matrix  can  be  inverted  by  first 
inverting  a  large  circulant  approximation  to  that  matrix,  followed  by 
inversion  of  a  much  smaller  Toplitz  matrix  [Ref.  1  ].  Hence  the  availa¬ 
bility  of  a  fast  optical  technique  for  inverting  large  circulant 
matrices,  together  with  a  digital  processor  that  inverts  a  small  Toplitz 
matrix,  allows  large  Toplitz  matrices  to  be  inverted,  hopefully  with 
greater  speed  than  afforded  by  an  all-digital  approach.  Once  a  method 
for  diagonalizing  circulant  matrices  optically  is  in  hand,  then  methods 
exist  for  inverting  such  matrices,  as  we  describe  later. 


There  is  an  additional  reason  for  interest  in  the  problem  of 
inverting  circulant  matrices  optically.  In  working  with  matrices, 
rather  than  pictorial  data,  one  is  dealing  with  numbers,  and  the  results 
of  the  processing  operation  are  likewise  numbers.  It  therefore  becomes 
rather  easy  to  assess  the  accuracy  of  the  operations  one  is  performing, 
and  to  discover  the  primary  sources  of  accuracy  limitation  in  the  sys¬ 
tem.  Thus  we  regard  this  project  as  providing  a  testbed  within  which 
the  sources  of  inaccuracy  of  coherent  optical  processors  can  be  studied 
in  a  more  precise  way  than  has  been  possible  before.  Knowing  the  chief 
reasons  for  inaccuracies  allows  one  to  focus  attention  on  these  sources 
of  error  and  to  explore  methods  for  reducing  their  effects.  Thus  infor¬ 
mation  discovered  in  this  project  may  have  wider  application  in  the 
field  of  coherent  optical  processing. 

Vith  this  information  as  background,  we  return  to  a  more  detailed 
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discussion  of  the  optical  problem.  The  task  of  diagonalizing  a  circulant 
matrix  using  coherent  optics  can  be  performed  if  the  normal  continuous 
„  two-dimensional  Fourier  transform  so  easily  performed  by  such  systems 

can  be  changed  to  a  discrete  Fabler  transform.  Much  of  our  work  on  this 

*  project  during  the  past  year  has  been  aimed  at  conversion  to  a  discrete 
Fourier  transform.  The  method  used  for  this  conversion  is  as  follows. 
The  matrix  to  be  diagonalized  is  entered  into  the  coherent  optical 
system  as  an  array  of  transmitting  cells  in  a  mask.  The  matrix  is 
repeated  at  least  3  by  3  times  in  the  horizontal  and  vertical 
directions,  causing  the  spectrum  to  form  a  series  of  discrete  spots. 
Each  of  these  spots  represents  a  different  complex  eigenvalue  of  the 
original  matrix.  Measurement  of  the  intensities  of  these  spots  by 
discrete  elements  of  a  detector  array  is  equivalent  to  measurement  of 
the  squared  magnitudes  of  the  eigenvalues  of  the  matrix.  If  the  full 
c<  mplex  values  are  desired,  then  interferometry  or  heterodyne  detection 
must  be  used  to  extract  both  amplitude  and  phase  information. 

A  major  component  of  our  effort  during  the  past  year  has  been  aimed 
at  an  analysis  of  this  method  for  performing  the  discrete  Fourier  trans¬ 
form.  The  analysis  examined  the  effects  of  the  repetitions  of  the 
matrix  at  the  input  plane,  the  finite  cell  sizes  used  to  represent 
t  matrix  elements,  and  the  finite  sizes  of  the  detector  elements  used  to 

measure  the  eigenvalues.  It  was  discovered  that  the  operation  performed 
by  the  two-dimensional  system  is  not  quite  that  indicated  in  Eq.  (3) 
above.  The  inverse  operation  associated  with  the  first  ¥  matrix  is 
missing  in  the  optical  realization,  and  as  a  consequence  the  locations 
of  the  spots  representing  the  eigenvalues  of  interest  are  not  quite 
those  expected  at  the  start.  Nonetheless,  the  desired  spots  are  present 
and  must  simply  be  detected  in  the  correct  region  of  the  output. 

Figure  1  below  shows  the  optical  setup  used  for  obtaining  eigen¬ 
values,  a  sketch  of  a  typical  matrix  transparency  introduced  into  the 
system  (for  a  4*4  case),  and  a  sketch  of  the  resulting  distribution  of 
light  intensity  in  the  focal  plane  of  the  lens.  Note  that  there  is  are 

* 

a  large  number  of  light  spots  in  the  output.  Those  representing  the 
eigenvalues  of  interest  are  enclosed  in  a  box. 


During  the  past  year  we  have  discovered  a  simple  technique  for 
entering  circulant  matrices  with  complex  elements  into  the  system. 
Matrices  representing  the  biased  real  part  and  the  biased  imaginary 
part  are  interlaced  with  a  diagonal  spatial  offset  of  one  quarter  of  a 
cell  separation  distance.  The  biases  do  not  affect  the  output  provided 
the  eigenvalues  are  detected  in  a  subdiagonal  that  does  not  contain  the 
aero- frequency  or  "D.C."  spot.  The  spatial  offset  of  the  real  and 
imaginary  part  matrices  results  in  the  complex  addition  of  their 
respective  eigenvalues  with  a  90  degree  phase  shift.  As  a  consequence 
the  eigenvalues  of  the  complex  matrix  are  detected.  The  technique  is 
illustrated  in  Figure  2. 

The  chief  accuracy  limitations  we  have  encountered  so  far  in 
measuring  eigenvalues  by  these  methods  have  arisen  from  inaccuracies  of 
the  optical  masks  representing  the  input  matrices.  These  masks  have 
been  written  on  a  DICOMED  plotter  at  the  NASA  Ames  Research  Center.  It 
was  necessary  for  us  first  to  calibrate  the  gray  levels  of  the  plotter. 
Following  this,  we  discovered  that  the  spot  sizes  written  by  the  plotter 
are  a  function  of  the  gray  level  being  written,  and  further  calibration 
had  to  be  done.  Finally  we  turned  to  having  the  masks  fabricated  in  our 
Integrated  Circuits  Laboratory.  Preliminary  results  with  one  such  mask 
show  that  moat  of  the  limitations  we  were  encountering  with  the 
transparencies  written  with  the  DICOMED  plotter  are  no  longer  present. 
In  our  latest  experiment,  using  only  a  3x3  matrix,  an  accuracy  of  1.535 
was  obtained  in  the  measured  eigenvalues.  This  accuracy  may  be  a 
function  of  matrix  size,  a  point  to  be  explored  in  the  future.  We  are 
also  obtaining  some  transparencies  written  with  an  electron  beam  litho¬ 
graphy  machine,  and  will  be  evaluating  them  for  accuracy. 

At  this  point  it  appears  that  our  major  problems  in  producing  masks 
that  accurately  represent  the  matrices  of  concern  have  been  overcome, 
and  that  future  efforts  can  be  more  directly  aimed  at  inaccuracies 
arising  within  the  coherent  optical  processing  system  itself. 

Attention  will  be  turning  later  this  year  away  from  the  measurement 


of  eigenvalues  towards  the  inversion  of  circulant  matrices.  In  order  to 
invert  such  matrices,  it  is  necessary  to  record  the  complex  eigenvalues 
in  the  focal  plane  holographically,  and  to  produce  a  transparency  with 
complex  amplitude  transmittance  proportional  to  the  reciprocal  of  the 
complex  fields  associated  with  the  original  eigenvalues.  A  second 
discrete  Fourier  transform  will  produce  an  output  matrix  that  is  the 
inverse  of  the  original  matrix.  The  necessary  nonlinearities  can  be 
achieved  with  photographic  film  using  the  techniques  that  have  been 
developed  in  the  past  for  coherent  optical  inverse  filtering.  We  intend 
to  do  our  work  with  such  methods,  although  our  interest  will  also  be 
focused  on  possible  techniques  for  performing  the  same  operations  in 
real  time  without  the  use  of  photographic  film.  One  such  technique  is 
the  subject  of  study  in  another  project  supported  by  this  grant  and 
described  in  the  section  that  follows. 

A  paper  describing  our  past  work  on  this  project  will  be  presented 
at  the  forthcoming  International  Optical  Computing  Conference. 

IV.  WAVEFRONT  DIVISION  BY  FOUR-WAVE  MIXING 

A  new  project  begun  during  the  past  year  is  aimed  at  using  four-wave 
mixing  to  divide  one  complex  wavefront  by  another  complex  wavefront. 
Past  work  by  others  has  indicated  that,  under  certain  conditions,  the 
strength  of  the  signal  generated  in  a  degenerate  four-wave  mixing 
experiment  is  proportional  to  the  contrast  of  the  fringes  between  the 
object  wave  and  one  of  the  pump  waves  [Ref.  2].  Workers  in  holography 
have  shown  that  under  such  conditions  holograms  can  be  used  to  perform  a 
division  between  a  complex  wavefront  incident  on  the  hologram  during 
reconstruction  and  the  complex  object  wavefront  stored  in  the  hologram 
[Ref.  3J.  The  idea  carries  over  directly  to  four-wave  mixing,  in  which 
the  wavefront  incident  in  the  probe  beam  will  be  divided  by  the  wave- 
front  carried  by  one  of  the  pump  beams.  The  chief  experimental 
difficulty  is  that,  with  currently  available  materials  (e.g.  BSO)  used 
in  the  way  necessary  to  achieve  division,  the  reflected  wave  of  interest 
is  extremely  weak.  This  problem  may  be  solved  with  the  availability  of 


crystals  with  higher  efficiency  and  even  gain,  hut  for  the  moment  the 
ideas  must  be  tested  with  BSO. 

A  major  effort  has  been  mounted  during  the  past  year  to  develop  the 
instrumentation  needed  to  measure  the  relevant  properties  of  BSO  at  the 
very  low  light  levels  anticipated.  This  work  has  involved  the  inter¬ 
facing  of  our  Reticon  array  to  a  microcompter  through  A/D  converters  in 
such  a  way  that  a  large  number  of  measurements  can  be  made,  the  results 
digitised,  and  the  measurements  averaged  to  increase  the  signal  to  noise 
ratio.  At  this  point  the  interfacing  is  done,  but  the  Reticon  array 
electronics  are  being  repaired  after  a  failure.  It  should  not  be  too 
long  until  we  can  take  the  measurements  needed  to  ascertain  the  dynamic 
range  over  which  wavefront  division  can  be  achieved. 

In  the  meantime  we  have  been  developing  ideas  regarding  possible 
applications  of  real-time  wavefront  division  to  practical  problems.  One 
fairly  obvious  application  is  to  inverse  filtering  or  image  deblurring, 
and  indeed  the  method  could  be  used  for  inversion  of  circulant  matrices, 
and  discussed  in  the  previous  section.  We  have  also  developed  some 
ideas  about  how  to  use  the  method  in  the  testing  of  integrated  circuit 
masks  for  defects,  but  it  is  a  bit  too  early  to  discuss  these  ideas  at 
this  point. 

As  part  of  our  effort  to  become  more  familiar  with  the  mathematics 
of  four-wave  mixing,  we  undertook  a  study  of  what  diffraction  efficien¬ 
cies  could  be  achieved  when  the  object  beam  contains  a  speckle  pattern, 
as  it  would  for  any  diffuse  object.  Relevant  work  in  holography  was 
found  [Ref.  4],  but  it  was  discovered  that  an  important  effect  was  left 
out  of  that  analysis  that  should  have  been  included.  This  effect  arises 
from  the  random  tilts  of  the  wavefront  (associated  with  the  speckle 
phenomenon)  and  the  resulting  random  tilts  of  the  fringes  of  the  inter¬ 
ference  patterns  in  the  crystal.  These  tilts  affect  the  diffraction 
efficiency  that  can  be  achieved.  It  waa  necessary  to  develop  further 
information  on  the  statistics  of  wavefront  tilt  in  q  speckle  pattern. 
Such  an  analysis  was  carried  out,  and  based  on  it,  a  paper  waa  aubmitted 
to  the  Journal  of  the  Optical  Society  of  America.  This  paper  has  now 
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been  accepted  and  will  appear  in  the  July  issue.  A  copy  is  appended  to 
this  report. 

The  major  emphasis  of  our  work  in  the  immediate  future  will  be  aimed 
at  collecting  the  necessary  experimental  data  that  will  tell  us  informa¬ 
tion  concerning  the  validity  of  the  ideas  about  wavefront  division. 
After  the  limitations  of  the  method  are  fully  understood,  attention  will 
be  devoted  to  studies  of  potential  applications  of  the  method. 


V.  STUDIES  0?  SPECKLE  SUPPRESSION  ALGORITHMS 


Another  new  effort  initiated  during  the  past  year  has  been  in  the 
area  of  speckle  suppression  in  coherently  formed  images.  Speckle  arises 
whenever  an  object  containing  surface  roughness  on  the  scale  of  the 
illuminating  wavelength  is  imaged.  Speckle  occurs  in  synthetic  aperture 
radar  images,  in  medical  -ultrasound  images,  and  milimeter  waves  images, 
for  example.  Any  method  that  would  suppress  this  disturbing  noise 
without  while  preserving  resolution  for  the  desired  object  information 
would  constitute  a  major  contribution. 

Many  methods  for  suppressing  speckle  have  been  studied  in  the  past. 
The  moat  obvious  method  is  simply  to  blur  the  intensity  distribution  in 
the  image  with  some  kind  of  linear  filter,  in  which  case  not  only  is 
speckle  reduced,  but  also  the  desired  portion  of  the  image  is  blurred 
with  a  consequent  loss  of  resolution.  The  key  question  to  be  addressed 
is  whether  there  exists  another  approach  to  speckle  suppression  that 
does  better  than  simple  blurring  in  terms  of  the  tradeoff  between 
speckle  noise  and  resolution.  Unfortunately,  many  recent  studies  of 
speckle  suppression  have  not  compared  the  results  obtained  with  those 
achievable  by  simple  blurring,  and  therefore  it  is  not  fully  possible  to 
judge  the  merits  of  these  techniques. 

In  beginning  our  work  in  this  area,  it  was  first  necessary  to 
develop  accurate  methods  for  simulating  the  speckle  phenomenon 
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digitally.  The  problem  is  far  more  difficult  than  might  be  imagined  at 
first  glance.  Simple  random  phase  coding  of  the  object,  followed  by 
low-pass  filtering  in  the  frequency  domain,  yields  as  serious  problem 
with  aliasing.  Our  approach  has  been  to  convolve  a  randomly  phase-coded 
object  with  an  accurate  repreaentation  of  the  amplitude  spread  function 
of  the  imaging  system,  without  resorting  to  discrete  Fourier  transforms 
at  all.  Excellent  results  have  been  obtained  for  one  dimensional 
objects,  and  we  are  now  extending  the  method  to  two  dimensions. 

The  speckle  suppression  methods  studied  (so  far  primarily  in  one 
dimension)  have  been  the  following: 

1.  Simple  blurring  of  the  intensity  image  by  a  rectangular  spread 
function. 

2.  Filtering  of  the  image  intensity  with  a  linear  Wiener  filter 
designed  for  multiplicative  noise. 

3.  Logarithmic  or  homomorphic  filtering,  in  which  the  logarithm  of  the 
image  intensity  ia  smoothed  with  a  rectangular  point-spread 
function,  following  which  exponentiation  of  the  image  is  performed. 

4.  Median  filtering,  in  which  the  median  of  the  image  intensity  values 
within  a  moving  rectangular  window  is  chosen  as  the  value  of  the 
filtered  image  at  each  point.  Several  iterative  passes  through  such 
filter  ultimately  yield  a  so-called  "root"  image  [Ref.  5],  which 
likewise  has  been  studied  as  a  possible  restoration  of  the  speckled 
image. 

5.  Maximum  entropy  processing,  an  approach  not  studied  before  in 
speckle  suppression,  and  adapted  to  this  problem  by  us.  More 
explanation  of  the  method  is  presented  below. 

In  all  cases  we  used  the  mean-squared  error  between  the  filtered 
image  and  an  "ideal"  image  as  the  criterion  of  goodness.  For  the 
"ideal"  image  we  considered  both  the  intensity  distribution  in  the 


coherently  formed  image  (obtained  from  a  8yatem  with  the  same  limiting 
pupil)  and  the  intensity  distribution  of  the  incoherently  formed  image 
(again  through  the  same  system).  The  results  indicated  that  all 
restoration  techniques  produce  images  that  are  closer  to  the 
incoherently  formed  image  than  to  the  speckle-free  coherently  formed 
image. 

The  maximum  entropy  approach  to  speckle  suppression  is  based  on  the 
following  ideas.  The  image  intensity  has  a  Fourier  transform  that  is 
identically  equal  to  the  deterministic  autocorrelation  function  of  the 
fields  incident  on  the  pupil  of  the  imaging  system.  Thus  the  problem  of 
estimating  the  image  intensity  in  the  presence  of  speckle  can  be 
regarded  as  a  kind  of  spectral  estimation  problem,  in  which  we  are 
trying  to  estimate  the  power  spectral  density  of  a  random  process  from 
one  sample  function  of  that  process.  If  we  were  able  to  exactly  estimate 
the  power  spectral  density  of  the  underlying  random  process,  we  would  be 
able  to  Fourier  transform  that  result  and  obtain  a  speckle-free  image 
with  full  the  resolution  of  the  imaging  system. 

At  low  spatial  frequencies  (small  separations  in  the  pupil  auto¬ 
correlation  process),  a  great  deal  of  spatial  averaging  is  performed, 
and  the  spectral  density  obtained  from  a  single  sample  function  (i.e.  a 
single  speckled  image)  is  very  close  to  the  value  that  would  be  obtained 
from  an  average  over  an  ensemble  of  images.  However,  at  high  spatial 
frequencies  relatively  little  spatial  averaging  occurs,  and  the  power 
spectrum  associated  with  the  single  speckled  image  departs  significantly 
from  that  associated  with  an  ensemble  average.  The  maximum  entropy 
approach  simply  throws  away  the  high  spatial  frequencies  in  the  spec¬ 
trum,  retaining  the  low  spatial  frequencies,  and  attempts  to  estimate 
the  high-frequency  portion  of  the  spectrum  from  the  low  frequency  por¬ 
tion  using  the  maximum  entropy  method.  Details  about  this  method  will 
be  presented  in  a  future  publication;  for  the  present  the  above  descrip¬ 
tion  should  suffice. 

Our  results  to  date  indicate  that  the  best  restorations  (best 
meaning  least  mean-squared  error)  are  obtained  with  a  linear  Viener 


filter  based  on  the  multiplicative  model  of  speckle*  Second  in  perfor¬ 
mance  is  the  maximum  entropy  method.  A  very  close  third  is  simple 
linear  averaging  with  a  rectangular  window.  Rather  far  behind  in 
performance  are  both  homomorphic  filtering  and  median  filtering. 

Criticisms  that  can  be  leveled  at  these  results  are  two-fold. 
First,  the  results  to  date  are  one-dimensional,  and  the  ordering  of 
performance  could  conceivably  be  different  in  two  dimensions.  Second, 
only  a  limited  set  of  objects  have  been  studied,  and  the  results  may  be 
object-dependent.  We  believe  that  neither  of  the  criticisms  will  turn 
out  to  be  valid,  but  further  work  is  necessary  to  verify  this  belief. 

The  maximum  entropy  method  used  is  a  rather  crude  one,  and  we  intend 
to  test  other  more  sophisticated  versions  of  this  approach.  The  project 
is  ultimately  a  sufficiently  large  one  that  we  intend  to  seek  other 
support  to  carry  it  on.  Hopefully  by  Autumn  of  this  year  other  support 
will  be  in  hand.  However,  AFOSH  will  continue  to  receive  credit  for 
their  initial  support  of  the  work,  without  which  it  could  not  have 
begun. 


VI.  PAPERS  PUBLISHED  AND  MEETING  PRESENTATIONS 

We  summarise  in  this  section  the  various  publications  and  presenta¬ 
tions  made  during  the  past  year  reporting  on  work  supported  by  this 
grant. 

A.  Papers  Published 

J.W.  Goodman  and  Moon  Song,  "Performance  limitations  of  an  analog  method 
for  solving  simultaneous  linear  equations”,  APPLIED  OPTICS  21,  502- 
506  (1982). 

M.  Tur,  K-C  Chin,  and  J.  W.  Goodman,  "When  is  speckle  noise 
multiplicative?",  APPLIED  OPTICS  21,  1157-1159  (1982). 


J. V.  Goodman,  Architectural  development  of  optical  data  procesaing 

systems"  (invited),  J.  E.  E.  E.  AUSTRALIA  2,  139-149  (1982). 

B.  Papers  Accepted  for  Publication 

K.  Ochoa  and  J.V.  Goodman,  "Statistical  distribution  of  Ray-directions 

in  a  fully  developed  speckle  pattern"  (accepted  for  publication  in 
J.  Opt.  Soc.  Am.) 

C.  Oral  Presentations 


J.V.  Goodman,  "Architectural  development  of  optical  data  processing 
systems",  Conference  on  lasers  and  electro-optics,  Phoenix,  Arizona, 
May,  1982  (invited). 

J.V.  Goodman,  "Optical  data  processing,  past,  present  and  future", 
Conference  on  Electrooptics,  Tokyo,  Japan,  December  1982  (invited). 

Finally,  the  contributions  of  various  individuals  should  be 
explicitly  mentioned.  Prof.  L.  Hesselink  provided  advice  and  encourag- 
ment  for  various  graduate  students  involved  with  the  grant.  Dr.  Moshe 
Tur  provided  guidance  for  graduate  students  working  on  various  problems 
described  above.  Vork  on  optical  interconnections  has  been  performed 
primarily  by  Raymond  Kostuk.  Vork  on  diagonalication  and  inversion  of 
circulant  matrices  was  performed  by  Q.  Cao.  Vork  on  wavefront  division 
using  four-wave  mixing  was  performed  by  E.  Ochoa,  who  was  also  partially 
supported  by  an  IBM  doctoral  fellowship.  Vork  on  speckle  suppression  was 
performed  by  Rae-Hong  Park. 
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FIGURE  CAPTIOUS 


1.  (a)  Optical  setup  for  obtaining  eigenvalues  of  circulant  matrices. 

(b)  Sketch  of  the  matrix  transparency. 

(c)  Sketch  of  the  focal  plane  light  distribution. 

2.  Method  for  encoding  complex  matricea. 
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INPUT  OF  A  4X4  CERULANT  MATRIX 


APPENDIX  I 


Architectural  Development  of  Optical  Data  Processing  Systems 

J.  W.  GOODMAN 


SUMMARY  The  architectural  development  of  optical  data  processing  ayatems  ia  traced,  beginning  with  the 
foundation  eatabliahed  by  Abbe  and  ending  with  a  diacuaaion  of  new  trenda.  Some  apeculation  on  the  future 
of  the  field  ia  alao  included. 


1  INTRODUCTION 

The  hiatoricel  development  of  optical  data  proceaa- 
ing  architecturea  can  be  represented  in  the  form  of 
a  tree.  Like  a  tree,  the  field  haa  roota,  a  atrong 
truok,  major  branchea  that  aubdivide  into  amaller 
limbs  and  new  growth.  In  many  esses  Che  breachea 
intertwine.  In  thia  paper  we  deacribe  the  major 
atructurea  of  thia  tree  and  apeculate  about  future 
growth.  The  reader  may  wiah  to  refer  to  fig.  1, 
showing  the  structure  of  this  tree,  aa  we  proceed 
through  the  material  that  followa.  The  picture 
preeented  is  a  personal  and  subjective  one,  as  it 
must  be  in  any  evaluation  of  past  developments  and 
speculation  about  the  future  yet  to  come. 


2  THE  ROOTS  AND  THE  TRUNK 

The  lower  roots  of  the  optical  data  processing  tree 
ere  undoubtedly  the  work  of  E.  Abbe  in  the  late 
19th  century  (1).  Abbe  recognised  that  coherent 
optical  iawging  systems  accept  and  reject  various 
grating  structures  associated  with  an  object.  He 
thereby  introduced  the  notion  that  such  aysteam  act 
aa  spatial  filters,  although  thia  term  was  not 
used.  The  upper  roots  of  the  tree  ere  perhaps  beat 
attributed  to  F.  Zernike,  whose  invention  of  the 
phase  contrast  microscope  (2]  provides  the  moat 
remarkable  early  example  of  a  successful  attempt  to 
manipulate  the  spatial  spectrum  of  a  coherent 
object  for  useful  purposes. 

The  lower  trunk  of  the  tree  is  the  work  of  Marechal 
(3],  Taujiuchi  (4],  O'Neill  (3),  Lohmann  (6],  and 
others,  who  were  among  the  first  to  apply  coherent 
optics  to  optical  data  processing  problems,  includ¬ 
ing  image  dablurring  and  extraction  of  two  dimen¬ 
sional  signals  from  noise.  The  upper  trunk  of  the 
tree  ia  the  work  of  Cutrona,  Leith,  Palermo  and 
Porcello  (7)  at  the  University  of  Michigan,  who 
first  emphasised  the  generality  of  optical  pro¬ 
cessors  and  outlined  a  multitude  of  methods  for 
putting  such  systems  to  use  in  the  field  of  radar 
signal  processing. 

In  the  early  1960's  the  tree  forked  into  several 
branches.  He  will  trace  the  development  of  these 
branches  in  the  following  sections. 


Invited  paper  submitted,  to  The  Institution 
of  Radio  and  Electronics  Rngineers  Australia 
on  30  March  1902. 


Figure  1  The  tree  describing  architectures  for 

optical  data  processing. 

3  THE  SYNTHETIC  APERTURE  RADAR  BRANCH 

One  of  the  first  and  most  important  branches  of  the 
optical  data  processing  tree  is  devoted  entirely  to 
the  processing  of  signals  collected  by  synthetic 
aperture  radars. 

A a  shown  in  fig.  2  (a),  an  aircraft  carrying  a  radar 
transmitter  end  receiver  flies  a  straight-line  path 
over  terrain  of  interest.  Radar  pulses  a.  •  trans¬ 
mitted  norsMl  to  the  flight  path  and  the  returning 
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echo**  arc  mixed  with  a  highly  stable  local 
oscillator  in  the  radar  receiver.  A  point 
scatterer  on  the  ground  returns  a  series  of  echoes 
as  the  aircraft  flies  past.  Figure  2  (b)  shows  a 
typical  record  of  the  returned  signal  from  two 
point  scatterers  at  different  ranges,  as  recorded 
in  the  aircraft.  It  is  asauaed  that  the  radar 
transaits  short  pulses  (although  the  use  of  longer 
"chirp"  pulses  is  also  permitted).  Range  inform¬ 
ation,  plotted  vertically,  is  obtained  by  pulse 
echo  timing,  while  azimuth  information,  plotted 
horizontally,  records  the  complex  amplitudes  of  the 
fields  received  from  scatterers  at  each  range  from 
the  flight  path.  From  this  data  it  is  desired  to 
fora  an  iaage  of  the  terrain  reflectivity.  Note 
that  the  record  from  a  single  scatterer  takes  the 
fora  of  a  one-diaensional  zone  plate,  with  a  focal 
length  that  depends  on  the  distance  of  the  scatter 
from  the  flight  path. 
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Figure  3  Range  and  azimuth  iaage  locations. 
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Figure  4  An  early  optical  system  for  processig 

synthetic  aperture  radar  data.  Note 
the  output  slit. 
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An  improved  optical  system,  called  the  tilted  plane 
processor,  was  described  in  1972  by  Kozraa,  Leith  and 
Massey  [9].  As  shown  in  fig.  3,  an  anaaorphic 
telescope  system  is  used,  together  with  a  tilted 
input  plane  and  a  tilted  output  plane,  to  bring  the 
range  and  azimuth  images  planes  into  coincidence 
with  a  magnification  that  is  no  longer  range 
dependent.  The  need  for  the  output  slit  is  thus 
eliminated  and  significantly  better  processor 
performance  is  obtained. 

The  tilted  plane  processor  reswins  today  one  of  the 
most  sophisticated  processors  in  the  entire  field  of 
optical  data  processing.  It  has  demonstrated  a 
continuing  usefulness,  even  in  Che  face  of  improved 
capability  of  digital  processing  systems. 


Figure  2  (a)  Synthetic  aperture  radar  data 

collection,  (b)  Record  obtained  for 
two  scatterers. 


Mien  short  pulses  are  used,  the  iswge  in  the  range 
direction  lies  at  the  film  on  which  the  data  was 
recorded.  However ,  due  to  the  dependence  of  the 
«one-place  focal  lengths  on  range,  the  azimuth 
atructure  of  the  image  lies  in  a  tilted  plane 
behind  the  film,  as  shown  in  fig.  3.  The  purpose 
of  the  optical  processor  is  to  bring  the  range  and 
asimuth  images  into  coincidence,  as  can  be  accom¬ 
plished  by  swans  of  an  anamorphic  optical  system. 

An  early  version  of  such  a  processor  (8]  used  a 
conical  lens  at  the  film  plane  to  move  the  azimuth 
image  to  infinity,  a  cylindrical  lens  to  the  right 
of  the  film  to  wove  the  range  image  to  infinity  and 
a  spherical  lens  further  to  Che  right  to  bring  the 
two  infinitely  distant  images  back  to  coincidence 
(see  fig.  4).  However,  the  image  so-formed  has  a 
magnification  that  is  range  dependent  and  therefore 
an  output  slit  was  required.  As  the  input  data 
film  and  the  output  recording  film  arc  moved  in 
synchronism-,  an  image  of  the  reflectivity  of  the 
terrain  is  recorded. 
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Figure  5  The  tilted  plane  processor.  Note  the 
absence  of  an  output  slit. 


4  THE  OPTICAL  PATTERN  RECOGNITION  BRANCH 

A  second  msjor  branch  of  the  tree  grew  out  of  the 
potential  application  of  coherent  optical  systems  to 
pattern  recognition  problems.  The  origins  of  this 
branch  lie  in  the  work  of  Vender  Lugt  [10}  inl964, 
who  originated  and  demonstrated  the  first  highly 
practical  method  for  generating  complex  matched 


■JvA- 


Jot* maf  irfJPW^aat  mas  iXf. 


*  VJ  v*i  iiwit  A ;Ai  A  *  9  1'91  L.I 


•  »Uw*  It  4  MW  V  til<  *  Xbv 


filters  for  two-dimensional  object*.  The  optical 
filtering  system  i*  tha  standard  one  ahovm  in  fig. 
6,  in  which  input  ia  introduced  by  means  of  fila  or 
a  real -tine  light  valve,  the  first  lent  Fourier 
tranafona*  the  input  fields,  the  filter  placed  in 
the  focal  plane  modifies  the  Fourier  spectrua  of 
the  object  and  the  second  len*  inverse  transforms 
the  Modified  spectrusi  to  produce  a  filtered  output 
image.  Vender  Lugt's  important  contribution  was 
the  origination  of  the  idea  of  an  interferooetric- 
ally  recorded  (or  holographic)  Fourier-plane 
filter,  sdiich  could  control  the  complex  aaplitude 
transmittance  through  the  focal  plane  in  an 
extreaely  flexible  way.  In  particular,  a  filter 
can  be  generated  that  ha*  as  it*  aaplitude  trans¬ 
mittance  the  complex  conjugate  of  the  Fourier 
spectrum  of  an  object  from  which  the  filter  was 
made.  Thu*  if  a  certain  object  o(x,y)  haa  as  it* 
Fourier  spectrua  the  complex  function  0(fx,f_),  a 
filter  with  amplitude  transmittance 

t(x, y)  -  0*(fx,fy)  (1) 

can  be  constructed.  Such  a  filter  is  said  to  be 
"matched"  to  the  object  o(x,y);  it  will  produce  a 
bright  output  spot  of  light  at  any  position  where 
the  particular  object  0(x,y)  is  present  ia  the 
input  field.  The  Vender  Luge  filter  has  played  an 
important  role  in  the  development  of  optical 
matched  filtering  approaches  to  pattern 
recognition. 


Figure  6  Standard  coherent  optical  processing 

system. 


While  the  idea  of  the  interferometrically  generated 
filter  for  matched  filtering  applications  has  had 
an  extremely  large  intellectual  impact,  it  must  be 
said  that  after  nearly  20  years  of  research  the 
applications  of  such  filter*  are  disappointingly 
limited.  The  difficulties  lie  not  with  the  concept 
of  the  interfsrometricelly  generated  filter,  but 
with  the  use  of  a  matched  filter.  While  the 
matched  filter  response  is  not  affected  by  pure 
translation  of  the  position  of  the  object  at  the 
input  (tha  output  spot  simply  moves  with  the 
input),  such  filters  are  exceedingly  sensitive  to 
scale  sise  end  rotation  of  the  object  to  be 
recognised.  That  is,  if  a  matched  filter  is 
constructed  for  recognition  of  an  object  with  one 
particular  scale  sis*  end  rotational  orientation, 
the  filter  will  usually  produce  very  little 
response  to  that  same  object  presented  at  the  input 
with  a  different  scale  sit*  or  orientation. 

Attempts  to  overcome  this  undesired  sensitivity 
have  been  clever  but  not  very  successful.  V* 
mention  in  particular  th*  elegant  Mellin  filtering 
approach  of  Casasent  and  Psaltis  Ill].  Th*  input 
data  is  presented  to  the  processor  in  polar  co¬ 
ordinates,  rather  than  rectangular  coordinates.  In 
addition,  th*  axis  corresponding  to  radial  position 
is  intentionally  strected  by  a  logarithmic  spatial 
distortion.  As  a  consequence,  it  can  be  shown  that 


the  Fourier  transform  of  the  radial  structure  of  the 
image  is  identical  with  a  so-called  Mellin  transform 
and  the  magnitude  of  the  Mellin  transform  ia 
invariant  under  changes  of  scale  site  (magnification 
of  the  input).  Simultaneously  the  optical  system 
Fourier  transforms  the  angular  structure  of  the 
input.  With  suitable  car*  to  take  account  of  edge 
effects,  rotation  of  the  object  ia  equivalent  to 
translation  in  the  angular  direction.  Since  switched 
filters  do  respond  properly  to  translated  versions 
of  the  structures  to  which  they  are  matched,  the 
Mellin  matched  filter  will  respond  properly  to  an 
input  object  regardless  of  its  scale  sis*  or  it* 
angular  rotation.  Unfortunately  the  process  of 
desensitising  the  system  to  magnification  and  rotat¬ 
ion  has  introduced  a  sensitivity  chat  was  not 
present  originally,  namely  to  the  position  of  the 
object  in  the  input  field.  The  Mellin  matched 
filter  will  not  respond  properly  to  objects  that 
have  a  different  position  in  the  input  field  than 
they  had  when  the  matched  filter  was  constructed. 
Casasent  and  Psaltis  [12]  solved  this  problem  by 
entering  into  the  system  only  the  power  apactmm 
of  the  input  pattern,  a  quantity  that  is  independent 
of  the  position  of  that  pattern.  However,  the 
discarding  of  all  phase  information  in  the  amplitude 
spectrum  results  in  a  decrease  in  the  ability  to 
discriminate  between  different  patterns.  As  a 
consequence  the  problem  of  siauiltaneous  position-, 
size-,  and  rotation-invariant  optical  pattern 
recognition  has  not  yet  been  solved  in  a  fully 
satisfactory  way. 
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Figure  7  Diffraction  pattern  sampling  system. 


The  "pattern  recognition"  branch  has  a  sub-branch 
that  is  particularly  worthy  of  mention.  We  refer  to 
so-called  "diffraction  pattern  sampling"  systems 
which  base  pattern  recognition  decisions  on  the 
structure  of  the  optically  obtained  Fourier  trans¬ 
form  of  an  input  [13].  The  systehi  is  illustrated  in 
fig.  7.  The  input  is  presented  to  th*  system  as  a 
transparency.  The  coherent  optical  system  Fourier 
transforms  that  input,  displaying  an  intensity 
distribution  across  the  back  focal  plane  that  is  the 
power  spectna*  of  the  input.  A  detector  consisting, 
for  example,  of  combinations  of  rings  and  wedges,  is 
used  to  extract  a  reduced  amount  of  data  from  the 
power  spectrum.  Th*  powers  measured  by  these 
detector  elements  are  then  digitised  and  digital 
pattern  recognition  algorithm*  are  applied  to 
classify  the  input.  This  system  ha*  several  prop¬ 
erties  that  make  it  potentially  vary  useful.  First, 
ths  detector  in  the  focal  plane  performs  an  all- 
important  reduction  of  information,  presenting  th* 
digital  system  with  data  of  reduced  eosplexity  but 
still  retaining  sufficient  information  for  th*  par- 
ticuler  task  at  hand.  Second,  th*  system  combines 
optical  and  digital  computations  in  an  appropriate 
way,  with  each  system  doing  computations  for  which 
it  is  particularly  well  suited,  tn  spit*  of  th* 
many  nice  features  to  racomaend  this  approach  to 
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this  approach  to  pattern  recognition,  it  hae 
received  leea  u*e  than  it  deaervea.  Perhaps  the 
increasing  emphasis  on  robotics  and  coeputer  vision 
will  lead  to  a  renewal  of  interest  in  the 
diffraction-pattern  sampling  approach. 

S  THE  ACOUSTO-OPTIC  SIGNAL  PROCESSING  BRANCH 

A  third  major  branch  of  the  optical  data  processing 
tree  derives  its  strength  from  developments  in 
acousto-optic  device  technology.  The  branch  begins 
in  the  early  1960's  with  the  work  Rosenthal  [14], 
Slobodin  [15],  and  Arm,  Lambert  and  Veissman  [16] 
on  the  use  of  acousto-optic  devices  for  signal 
processing.  An  excellent  review  of  this  area  it 
found  in  reference  17  and  accompanying  papers  in 
the  same  issue. 


Figure  8 


(a)  Acousto-optic  cell  in  the  Raman- 
Nath  regime,  (b)  Acousto-optic  cell 

in  the  Bragg  regime. 


Before  beginning  a  discussion  of  acousto-optic 
signal  processing  architectures,  it  is  perhaps 
worthwhile  to  present  a  small  amount  of  background 
on  the  acousto-optic  devices  themselves.  Figure  8 
illustrates  acousto-optic  cells  in  two  different 
modes  of  operation,  in  both  cases,  an  electrical 
signal  v(t),  consisting  of  an  amplitude  and  phase 
modulated  carrier, 

v(t)  -  A(t)cos(2*v0t  ♦  (t)),  (2) 

is  applied  to  an  acousto-optic  transducer,  idiich 
launches  an  acoustic  wave  in  the  transparent  medium 
of  the  cell.  The  presence  of  the  acoustic  signal 
results  in  local  changes  of  the  refractive  index  of 
the  call,  with  the  result  that  an  optical  signal 
propagating  across  the  sound  beam  expariancaa 


spatially  varying  phase  modulation.  Figure  8  (a) 
shows  the  so-called  Raman-Nath  regime  of  operation, 
in  which  the  cell  thickness  and  spatial  period  of 
the  acoustic  signal  are  such  that  the  cell  acts  as  a 
"thin"  phase  grating,  generating  a  multitude  of 
diffraction  orders,  each  with  a  different  propa¬ 
gation  angle  and  a  different  optical  frequency;  only 
the  sero  and  first  order  beams  are  shown.  Figure  8 
(b)  shows  the  more  common  Bragg  regime,  in  which  the 
cell  acts  as  a  "thick"  phase  grating,  with  only  one 
or  two  diffraction  orders  of  significant  intens¬ 
ity  and  requiring  illumination  of  the  cell  at  the 
Bragg  angle  in  order  to  produce  a  strong  component 
of  diffracted  light.  Early  attention  to  devices 
operating  in  the  Raman-Nath  regime  soon  turned  to 
devices  operating  in  the  Bragg  regime,  idiere  wider 
bandwidths  can  be  achieved. 

This  acousto-optic  signal  processing  branch  has  had 
remarkable  growth  in  recent  years.  This  growth  has 
been  caused  partially  by  advances  in  the  technology 
of  acousto-optic  cells  and  partially  by  development 
of  new  architectures  in  which  such  cells  can  be 
used.  The  branch  can  be  said  to  have  forked  into 
two  major  sub-branches,  space-integrating  processors 
and  time-intergrating  processors.  Each  of  these 
sub-branches  has  a  complicated  branching  structure 
and  tht-  growth  is  so  dense  that  the  two  are  at  least 
partially  intertwined. 

Most  important  in  the  space-integrating  lineage  is 
the  Bragg-cell  spectrum  analyser,  shown  in  fig.  9, 
which  has  enjoyed  a  success  in  application  that  is 
beginning  to  rival  that  of  the  tilted-plene  pro¬ 
cessor  mentioned  earlier  (although  the  applications 
are  quite  different).  The  Fourier  transforming 
properties  of  coherent  optical  systems  areused  here 
to  perform  spectrum  analysis  of  temporal  signals, 
the  temporal  signals  being  introduced  into  the 
processor  as  spatial  signals  travelling  through  the 
input  acousto-optic  cell.  Coherent  optical  spectrum 
analysers  with  approximately  one  Chs  of  bandwidth 
and  time-bandwidth  products  of  approximately  1000 
are  now  in  use  in  a  number  of  laboratories  and  the 
number  of  lower  performance  systems  is  even  more 
plentiful.  Integrated  optic  versions  of  such 
spectrum  analysers  are  also  being  constructed,  based 
on  surface-acoustic-wave  input  of  RF  signals. 
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Figure  9 


INPUT 


Bragg  spectrum  anaelyeer. 


A  second  architecture  in  this  seme  space-integrating 
lineage  is  the  space-integrating  correlator,  in 
which  one  or  even  two  Bragg  cells  nay  be  used  to 
perform  correlation  of  two  signals.  A  typical 
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configuration  it  thown  in  fig.  10.  Tht  r. f.  aignal 
V|(t)t  having  complex  envelope 

Ait)  •  |  Ait)  |  exp(j0(t)|,  i3) 

where  Ait)  and  0(t)  are  the  amplitude  and  phaae 
■odulationa  of  vi(t),  it  to  be  correlated  with  a 
tecond  tignal  vjix)  (complex  envelope  Bix))  whichit 
ttored  at  a  fixed  reference  function  on  a  mask. 
Oting  the  optical  system  thown  in  the  figure,  the 
current  generated  by  a  point'  detector  located  oo 
the  optical  exit  in  the  output  plane  taket  the  fora 

i*it)  -  |  [  A* (x-Vt)  Bix)  rect{x/W)  dx  |  2  (4) 


where  V  it  the  velocity  of  propagation  of  the 
acouatic  wave  in  the  input  cell,  W  it  the  length  of 
the  illuainated  region  of  the  acouatic  cell  and 
reetix)  it  defined  to  be  unity  for  and  tero 

otherwite.  At  the  naae  iapliee,  the  correlation 
integral  ia  perforaed  over  apace,  while  different 
valuet  of  delay  are  realiaed  in  tine. 
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by  the  length  of  the  tine  integration,  which  in  turn 
it  liaited  by  the  dynaaic  range  of  the  detector. 

The  apace-bandwidth  product  of  the  procettor  liaitt 
the  range  of  achievable  delaya  for  which  the  crosa- 
correlation  function  can  be  awaaured. 


Figure  11  Tiae-integrating  correlator. 


An  alternate  fora  of  the  tiae-integrating  correl¬ 
ator,  at  introduced  by  Kellaan  (20]  and  Turpin  (21] 
ia  thown  in  fig.  12.  In  thia  cate  one  of  the 
aignala  ia  introduced  aa  a  tine  modulation  of  the 
intenaity  of  an  LED  or  later  diode  aource,  while  the 
tecond  aignal  ia  introduced  via  an  acouato-optic 
cell.  The  apatial  filter  thown  in  the  focal  plane 
of  the  tecond  lent  blocka  all  but  the  plua  firat  and 
aero  diffraction  ordera  and  ahifta  the  phaae  of  the 
aero  order  by  90  degreea  with  reapect  to  the  firat 
order. 


Figure  10  Space-integrating  correlator. 


The  tub-branch  repreaenting  time-integrating 
corrclatora  begin*  in  1972  with  the  patent  of  R.W. 
Hoot gome ry  (18|.  A  more  widely  known  developmentof 
these  ideas  ia  found  in  the  work  of  Sprague  and 
Koliopoulos  (19]  in  1976.  Thie  architecture,  which 
ia  illustrated  in  fig.  11,  uses  integration  in  time 
for  performing  the  correlation  integral  itself, 
•diile  using  apace  to  represent  various  values  of 
delay.  One  of  the  aignala  to  be  correlated,  vi(t), 
ia  input  a*  a  traveling  acouatic  wave  in  the  first 
cell,  while  the  second  aignal  ia  introduced  at  a 
counterpropagating  wave  in  the  second  cell.  The 
Bragg  effect  results  in  only  a  tero  order  and  a 
negative  first  order  being  transmitted  by  the  firat 
cell,  idiile  each  of  these  orders  ia  split  into  a 
aero  order  and  a  negative  first  order  by  the 
second  cell.  The  tero-order  beams  are  both 
blocked.  The  first-order  beams,  tdiich  are  at  the 
•aw  optical  frequency,  are  both  paesed  and  are 
ellowed  to  interfere  on  e  tint  integrating  detector 
array  at  the  output.  The  apatial  distribution  of 
time  integrated  intenaity  takes  the  fora 

Rg(x)  -  bias  term* 

♦  Be  /  [  A*(t-x/V)B<t*x/V)  dt)  (5) 


where  T  represent*  the  total  integration  time 
(liaited  in  practice  by  bias  build-up).  Mote  that 
for  each  different  value  of  x,  the  integrated 
latensity  ia  determined  by  the  value  of  the  cross- 
correletion  function,  evaluated  et  time  delay  2x/V. 
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Figure  12  Time-integrating  correlator  with  a 
modulated  source. 


A  close  relative  of  the  time-integrating  correlator 
is  the  triple  product  processor  of  Kellman  (22], 
which  allow*  the  more  general  class  of  two- 
dimensional  operation* 

Kd(x,y)  ■  bias  term* 

♦  j  vj(t)  vjft-x/V)  vj(t-y/V)dt  (6) 

T 

to  be  performed. 

The  acouato-optic  aignal  processing  branch  is  an 
axtremely  healthy  and  vigorous  one,  with  continued 
growth  at  the  tip*  of  the  branches.  There  has 
recently  been  a  tendancy  for  the  space-integrating 
and  time- integrating  branches  to  grow  back  together, 
with  processor*  being  proposed  that  use  both  types 
of  integration  simultaneously. 

6  DISCRETE  OPTICAL  PROCESSORS 


The  significance  of  the  time-integration  architect¬ 
ure  is  that  the  time-bandwidth  product  achieved  is 
no  longer  constrained  by  the  apace-bandwidth 
produet  of  the  processor  but  rather  ia  determined 


A  relatively  young  vigorous  branch  of  the  architect¬ 
ural  tree  represent*  activities  in  the  development 
of  discrete  optical  processors.  Such  processor*  can 
perform  the  discrete  analoge  of  any  continuous-time 

—  _  ■  r 


1 intar  operation*.  The  general  type  of  operation 
perforate!  ia  a  matrix-vector  product  of  the  fora 

y  -  Ax 

where  x  i»  a  length  M  input  coluan  vector,  A  ia  an 
NxH  matrix  and  y  ia  a  length  H  output  column 
vector.  It  ahould  be  eaphaaiaed  that,  while  the 
ayateaa  of  concern  are  diacrete  (i.e.,  baaed  on 
diacrete  mathematics),  they  are  none-the-leaa 
analog  aystema  and  subject  to  the  uaual  limitation* 
on  accuracy  and  dynamic  range  typical  of  auch 
ayateaa. 

Early  veraion*  of  optical  matrix-vector  multiplier* 
ware  baaed  on  the  uae  of  coherent  light  (23-5). 
Incoherent  implementation*  were  alao  deviaed 
(26,27) . 

However,  a  moat  aignificant  advance  in  thia  area 
came  from  the  work  of  Bocker  (28),  Bromley  (29)  and 
Monahan  et  at  (30),  who  deviaed  a  method  for 
matrix-vector  multiplication  uaing  a  tingle 
Light-emitting  diode  (LED)  aource  and  a  two  dinens- 
ional  charge-coupled-device  (CCD)  detector  array, 
aa  illustrated  in  fig.  13.  Thia  development  i* 
represented  by  an  entire  aub-branch  in  our  tree. 

In  explaining  the  baaic  operation  of  such  a  aystem, 
we  assume  for  Che  moment  that  all  elements  of  the 
input  vector  x  and  Che  matrix  A  are  non-negative 
and  real,  ao  that  they  can  be  represented  by  light 
intenaitie*  or  intensity  transmittance*.  The 
eleawnt*  of  the  input  vector  x  are  entered  into  the 
system  in  time  sequence  aa  pulse*  with  intensities 
proportional  to  the  desired  values.  The  light  from 
each  pulse  diverges  and  fall*  on  the  entire  matrix 
mask,  tdiich  contains  an  array  of  calls  with  intens¬ 
ity  transmittance*  proportional  to  the  element*  of 
the  desired  matrix.  The  light  transmitted  by  the 
suak  fall*  on  the  two-dimensional  CCD  detector. 

The  clocking  of  the  CCD  device  is  auch  that  charges 
accumulated  in  a  given  cell  are  tranafered  one 
column  to  the  right  (i.e.,  out  of  the  paper  in  the 
figure)  along  a  fixed  row  before  the  occurrence  of 
the  next  optical  pulae  and  in  this  way  the  contri¬ 
butions  of  all  pulses  are  accumulated,  row  by  row. 
For  example,  if  the  first  pulse  represents  the 
vector  element  xj,  and  hk,j  represent*  the  element 
of  the  matrix  mask  in  the  kth  row  and  first  column, 
the  cherge  deposited  in  the  kth  detector  element  of 
the  first  colunn  is  proportional  to 

Qk(l>  "  >>k»lxl* 

This  charge  is  clocked  one  column  to  the  right  and 
the  second  pulse,  with  intensity  proportional  to 
the  vector  element  *2  is  emitted.  In  the  second 
column  and  the  kth  row,  the  total  accumulated 
charge  becomes 

Qk<J>  -  h*,^  ♦  hk,2K2. 


After  M  pulses  and  M  charge  transfers,  the  total 
charge  accumulated  in  the  kth  element  of  the  last 
coluan  of  the  detector  is 


X-l 

QfcfM)  •  Z  hk_Xp 

p-0 

which  is  precisely  the  kth  element  of  the  desired 
output  vector  y.  In  principle,  all  N  elements  of 
the  output  vector  are  available  in  parallel  but 
only  after  M  cycle*  of  the  clock. 
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Figure  13  Serial  incoherent  matrix-vector 
multiplier. 


A  system  of  the  type  described  above  is  capable  of 
multiplying  a  matrix  of  approximately  500  x  500 
elements  times  a  500-length  vector  of  non-negative 
and  real  elements  with  clock  rates  of  approximately 
lOMHs. 

i 

Incoherent  system*  of  this  kind  are  not  restricted  j 
to  performing  operations  involving  only  non-negative 
and  real  quantities.  It  is  possible  to  code  any 
bipolar  real  quantity  and  indeed  any  complex-valued 
quantity  in  term*  of  two  nonnegative  and  real  quan¬ 
tities,  although  in  the  latter  case  certain  bias 
terms  must  be  subtracted  from  the  computed  result 
(31).  He  do  not  dwell  on  the  details  here,  it  being 
sufficient  to  say  that  fully  coeplex  operations  can 
be  performed  at  the  price  of  doubling  the  sire  of 
the  awtrix  and  the  length  of  the  vectors. 

A  second  sub-branch  of  the  branch  of  discrete 
optical  processors  is  represented  by  a  closely 
related  but  faster  incoherent  discrete  system  intro¬ 
duced  by  Goodman,  Dias  and  Woody  (32)  in  1978.  In 
this  case,  as  illustrated  in  fig.  14,  the  elements 
of  the  input  vector  are  introduced  in  parallel  on  ad 
array  of  LED's  and  the  elements  of  the  output  vectod 
are  detected  in  parallel  on  an  array  of  independent 
detector  elements.  The  first  box  labeled  "optics" 
spreads  the  light  from  each  LED  into  a  vertical 
column,  tdiich  illuminate*  a  single  coluam  of  the 
aMtrix  mask.  The  second  box  labeled  "optics"  placml 
on  aach  photodetector  a  light  intensity  proportional 
to  the  sum  of  the  light  intensities  transmitted  by 
an  entire  row  of  the  matrix  auk.  These  two  operat 
ions  result  in  detected  signals  proportional  to  thm 
desired  eleawnts  of  the  output  vector  y.  Versions 
of  this  system  based  on  discrete  optical  componentt 
(i.e.,  anamorphic  systems  of  lenaaa)  and  multi-mode 
optical  waveguides  have  been  built  (33).  The  chief 
disadvantages  of  the  system  are  the  relatively  shoe 
lengths  of  the  vectors  Chat  can  be  accommodated  and 
the  absence  of  any  data  compression  or  information 
reduction  in  the  operations  performed.  The  latter  I 
fact  result*  in  exceedingly  high  rates  with  which  1 
data  pours  out  of  the  system  (potentially  as  high  I 


as  10  GHx  or  more  for  a  100-channel  system) 


Figure  14  Parallel  incoherent  matrix-vector 
multiplier. 


An  exceedingly  intereacing  development,  const  it  ut- 
ing  an  extenaion  of  the  eub-branch  diacuaaed  above, 
ia  the  uae  of  the  parallel  incoherent  matrix-vector 
multiplier  in  an  iterative  a»de  aa  a  meana  for 
aolving  aimultaneoua  aeta  of  linear  equationa 
(34).  In  thia  caae  ue  ara  given  a  aet  of  simultan- 
eoua  linear  equationa  of  the  form 

Hx  •  c 

where  M  ia  a  known  matrix  and  c  ia  a  known  vector 
and  our  goal  ia  to  determine  the  elementa  of  the 
unknown  vector  x.  A  feedback  arrangement  aa 
illustrated  in  fig.  IS  can  be  ahown  to  reault 
ideally,  after  a  aufficiant  number  of  iterations, 
in  convergence  of  the  output  vector  to  the  deaired 
solution  x.  This  innovation  ia  significant  for 
optical  data  processing  for  the  following  reason. 
Nearly  all  previous  uses  of  optical  processing  have 
been  for  performing  operations  such  as  convolut¬ 
ions,  correlations  and  Fourier  transforms,  where 
the  number  of  elementary  operations  required  is  of 
order  N2  (N  being  the  number  of  degrees  of  freedom 
of  the  input)  but  for  which  fast  algorithms  exist 
reducing  the  required  number  of  operations  for  a 
digital  implementation  to  order  RlogN.  The  problem 
of  solving  a  set  of  N  simultaneous  linear  equations 
for  the  values  of  N  unknowns  is  fundamentally  an 
operation  requiring  order  operations  and,  ia  the 
general  case,  no  fast  algorithms  exist.  Any  appli¬ 
cation  of  optical  processing  to  the  domain  requir¬ 
ing  N2  operations  has  potentially  great  rewards, 
for  in  such  cases  the  "crutch"  of  fast  algorithms 
is  no  longer  available  to  the  digital  hardware  that 
competes  with  the  optical  processor  and  the 
inherent  speed  and  parallelism  of  the  optical 
approach  has  its  greatest  significance. 


eigenvectors  of  suttrices  have  been  proposed  [3S,  36] 
and  no  doubt  other  clever  applications  of  this  type 
of  processor  will  be  found  in  the  future.  However, 
such  processors  have  not  yet  been  fully  analysed 
from  the  point  of  view  of  accuracy  and  numerical 
stability.  The  chief  disadvantage  of  optical  pro¬ 
cessors  in  these  applications  is  clearly  the  limited 
accuracy  associated  with  any  analog  approach  to 
computation.  This  disadvantage  will  undoubtedly 
limit  the  kinds  of  matrices  that  can  be  successfully 
dealt  with  in  this  fashion  (37). 

The  last  sub-branch  from  the  discrete  optical  branch 
is  one  we  will  call  "systolic  processors".  This 
sub-branch  is  exceedingly  young  and  isasature  but 
also  exceedingly  interesting  as  a  novel  approach  to 
optical  cos^utation.  The  only  published  works  on 
optical  systolic  processors  at  the  time  of  this 
writing  are  those  of  Caulfield,  Rhodes,  Foster  and 
Horvits  [38]  and  Casasent  [39]. 
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(a)  Basic  building  block  of  a  systolic 
processor,  (b)  Thrse  processors 
interconnected  in  a  systolic  array. 


Figure  13  Incoherent  matrix-vector  multiplier 
with  feedback. 


The  area  of  iterative'  optical  processing  remains  an 
active  one,  with  considerable  potential  for  useful 
contributions.  ‘Recently  optical  implementations  of 
iterative  methods  for  finding  the  eigenvalues  and 


The  systolic  processor  architecture  is  the  invention 
of  H.  T.  Rung  and  C.  S.  Leiserson  [40]  and  has  been 
considered  primari 1 y  as  an  architecture  for  VLSI 
implementation.  There  are  in  fact  many  different 
architectures  for  systolic  processors,  depending  on 
the  operation  to  be  performed.  He  consider  here 
only  the  simplest  of  these,  namely  one  designed  to 
perform  the  aMtrix-vector  product  discussed 
earlier.  The  basic  building  block  of  the  systolic 
processor  is  illustrated  in  fig.  16(a).  An  input  x, 
represanting  one  elesMnt  of  the  input  vector, 
arrives  from  the  left.  Simultaneously  a  weighting 
coefficient  a,  representing  an  element  of  the  matrix 
A,  arrives  froo  above.  The  processor,  represented 
by  the  box,  receives  the  input  x  from  the  left  and 
passes  it  unchanged  to  the  right  along  the  upper 


(input)  and  cha  CCD  device  (output),  with  the 
coupling  supplied  by  the  strobed  LCD  sources. 


input  line;  simultaneously  it  accepts  the  value  y 
coning  Iron  the  right  on  the  output  line  and 
transforms  it  into  the  new  value  y+ax,  which 
continues  to  pass  to  the  left. 

figure  16(b)  shows  three  basic  processors  connected 
together  in  a  structure  that  will  Multiply  a  2x2 
matrix  A  times  a  length  2  input  vector  x,  producing 
a  length  2  output  vector  y.  He  can  regard  this 
structure  as  being  a  pair  of  coupled  delay  lines  in 
which  input  and  output  signals  counter-propagate 
and  with  coupling  coefficients  that  depend  on  the 
elements  of  the  aatrix  A.  The  elements  of  the 
matrix  are  displayed  above  the  processors  and  are 
assumed  to  arrive  with  a  timing  represented  by 
their  vertical  distance  above  the  processors.  Any 
time  interval  represented  by  a  dot  in  this  diagram 
contains  a  coupling  coefficient  that  is  immaterial, 
since  it  will  not  affect  the  output  signals  of 
interest.  The  input  deta  is  entered  from  the  left 
as  a  series  of  values,  one  every  two  cytles.  The 
first  input  vector  element,  xj,  passes  through  the 
first  processor,  generating  an  output  propagating 
to  the  left  on  the  output  line.  This  value  is  not 
of  interest  and  is  ignored.  The  input  clement  then 
passes  to  the  second  processor,  arriving  coincident 
with  the  arrival  of  coupling  coefficient  a]|.  The 
result  is  an  output  value  propagating  to  the 

left  on  the  output  line.  This  value  arrives  at  the 
processor  on  the  far  left  at  the  end  of  one  cycle 
period.  Meanwhile,  a  second  input  clement  xj  ha* 
been  launched  on  the  input  line  and  arrives  at  the 
first  processor  simultaneously  with  the  arrival  of 
the  above  partial  output.  The  coupling  coefficient 
at  Che  first  processor  is  a^.  resulting  in  a 
signal  at  the  final  output  on  the  left  of 

n  “  *11*1  ♦  *12x2 

In  a  similar  way,  the  various  coupled  inputs  yield 
S  second  output  signal  of  the  form 


y2  -  aj^j  ♦  *22*2  • 

i.e. ,  the  second  element  of  the  output  vector  y. 
The  concept  can  clearly  be  extended  to  larger 
matrices  and  vectors. 
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Figure  17  Optical  systolic  matrix-vector 
aniltiplier. 


An  optical  analog  version  of  the  systolic  matrix- 
vector  multiplier  has  been  proposed  in  reference  38 
and  ia  illustrated  in  fig.  17,  again  for  the  simple 
case  of  a  2x2  matrix.  The  elements  of  the  input 
vector  are  entered  as  propogating  pulses  of  short 
duration  in  an  acousto-optic  cell.  The  matrix 
elements  are  input  in  proper  time  sequence  as 
intensity  modulations  of  LCD's  in  a  parallel 
•Cray.  The  output  elements  accumulate  as  moving 
charge  packets  on  a  clocked  CCD  device,  fed  by 
aignals  from  a  parallel  array  of  detectors.  The 
two  delay  lines  involved  are  the  acousto-optic  cell 


Another  version  of  an  optical  systolic  processor  is 
under  construction  at  Stanford  Diversity  (41].  The 
processor,  which  we  call  a  fibre -op  tic  scattering 
processor,  is  illustrated  in  fig.  18.  It  consists  of 
two  single-mode  optical  fibres  weakly  coupled  by  a 
series  of  fibre-optic  couplers.  Ideally  the 
couplers  should  be  changeable  with  time  but,  as 
described  shortly,  useful  signal  processing  operat¬ 
ions  can  be  performed  even  with  fixed  couplers.  The 
elements  of  the  input  vector  are  entered  as  pulses 
of  varying  intensities  on  the  input  fibre,  while  the 
output  elements  accumulate  as  they  propagate  along 
the  output  fibre  and  emerge  in  time  aequence.  If 
the  couplers  have  fixed  coupling  coefficients, 
rather  than  being  time  changeable,  the  structure  of 
the  aMtrix  A  is  effectively  constrained  to  be  of 
Toeplits  form,  i.e.,  the  elements  along  any  one 
subdiagonal  are  all  identical.  Such  forms  occur 
when  the  sMtrix-vector  operation  represents  a 
convolution  or  a  correlation,  rather  than  a  more 
general  time-variant  operation. 
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Figure  18  The  fibre  optic  scattering  processor. 


It  is  difficult  to  avoid  the  feeling  that  these  new 
architectures  are  but  the  beginning  of  a  branch  that 
will  soon  be  developing  a  complicated  network  of 
offshoots,  most  of  which  are  impossible  to  fully 
visualise  or  iswgine  at  this  particular  point 
intime. 

7  SOME  SPECULATION  ABOUT  THE  fUTURE 

He  have  surveyed  the  past  and  present  architectures 
of  optical  data  processing  systems.  Perhaps  in 
closing  we  will  be  permitted  to  speculate  a  bit 
about  the  future. 

A  reasonable  starting  point  for  this  speculation  is 
considaration  of  tha  weak  and  strong  points  of 
optical  processors.  A  chief  weakness  of  such  data 
processing  systems  is,  in  our  view,  their  inherent 
lack  of  accuracy,  arising  from  the  fact  that  they 
are  purely  analog  in  nature.  A  further  weakness  is 
the  limited  number  of  different  types  of  basic 
operations  such  systems  can  perform,  generally 
additions,  subtractions,  and  multiplications.  The 
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nd««l«|  features  or  tnt  ©ptleai  pVoeeiiiotTStM 
•pood  with  which  computations  con  bo  performed  in  • 
highly  parallel  fashion. 

By  comparison,  digital  electronic  proceeeore  can 
achieve  accuracioe  liaited  only  by  the  number  of 
bice  retained  in  the  coaputotion.  In  addition, 
their  repertoire  of  type*  of  basic  operations  is 
nearly  unliaited. 

The  epeed  advantage  associated  with  optics  needs 
closer  axaalnation.  The  input  and  output  devices 
used  with  such  systeas  are  usually  electronic  in 
nature  and  therefore  there  ie  generally  an  elec¬ 
tronic  bottleneck,  liaiting  the  rate  at  tdiich  data 
can  enter  and  exit  the  aystea.  These  electronic 
bottlenecks  are  partially  coapensated  for  by  the 
parallalisa  with  which  the  analog  coaputations  are 
perforaod  within  the  proceseor.  The  root  of  this 
inherent  parallelisa  lies  in  Che  relative  ease  with 
trttich  optical  signals  can  be  used  to  provide 
aultiple  independent  paths  between  different 
coaputational  coaponents  of  the  processor  (e.g., 
transparencies  and  lenses).  Thus  the  strength  of 
optics  lies  in  its  ability  to  provide  intevaameat- 
icno  of  an  extreaely  coaplex  type  between  different 
coaputational  aleamnts,  while  the  strength  of 
digital  electronics  lies  in  the  ability  of  its 
coaputational  eleaenta  to  provide  flexible  operat¬ 
ions  with  high  accuracy. 

The  question  now  quite  naturally  arises  as  to 
whether  it  night  be  possible  to  combine  the  good 
features  of  optics  and  electronics  in  a  single  pro- 
ceesor.  We  suggest  that  the  answer  in  the  future 
will  be  affiraative.  We  envision  a  processor  that 
is  basically  an  electronic  chip  or  a  aeries  of 
chips,  with  digital  electronic  coaputational 
eleaenta  but  with  interconnections  between  elements 
and/or  between  chips  provided  by  optics.  While  the 
idea  of  using  optical  waveguides  on  a  chip  may  have 
asrit  in  some  situations,  it  is  our  view  that  there 
is  too  little  difference  between  aetallic  inter- 
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connections  and  waveguides  for  this  change  to  be  of 
paramount  significance.  Rather  we  would  speculate 
that  imaging  interconnections  are  worthy  of 
consideration.  To  illustrate  what  we  mean  by  this 
term,  consider  the  geometry  shown  in  fig.  19(a),  in 
which  two  chips  are  communicating  by  aeans  of 
isMging  interconnections.  At  various  points  on  the 
chip  on  the  left,  coaputations  have  been  completed 
end  are  ready  for  transfer  to  various  points  on  the 
chip  on  the  right.  The  chip  on  the  left  is  assumed 
to  contain  several  tiny  sources  of  light  or,  in  the 
case  of  light  supplied  froa  off  the  chip,  several 
tiny  andulators.  The  chip  on  the  left  ie  assumed 
to  contain  several  detectore  at  different  locat¬ 
ions.  Bach  source  (or  modulator)  on  the  left  is 
then  imaged  onto  an  appropriate  set  of  detectors  on 
the  right.  The  imaging  element  ie  certainly  more 
complicated  then  a  simple  lens.  It  might  be  a 
holographic  optical  element. 

Of  the  many  problems  with  the  above  idea  that  eoae 
to  adnd,  chief  among  these  is  perhaps  the  fact  that 
we  do  not  know  how  to  integrate  optical  sources  on 
silicon.  As  indicated,  however,  it  is  not  actually 
necessary  to  have  sources  on  the  chip,  only 
siodulators  that  can  modify  externally  supplied 
light.  Nonetheless,  these  modulators  aust  be 
exceedingly  small  in  sise  and  capable  of  high-speed 
operation  with  very  small  power  dissipation. 

The  ideas  presented  above  can  be  extended  to  the 
case  of  a  single -chip .as  shown  in  fig.  19(b).  In 
this  case  the  holographic  imaging  element  is 
reflective  and  resides  above  the  chip. 


Figure  19  (a)  Imaging  interconnections  between 

two  chips  by  means  of  a  holographic 
optical  element,  (b)  Imaging  inter¬ 
connections  within  a  single  chip. 


While  the  possibility  of  having  sources  or  modulat¬ 
ors  on  the  chip  may  seem  remote  at  this  point  in 
time,  there  is  a  subset  of  problems  for  which  they 
are  not  needed  in  integrated  form.  This  subset 
consists  of  data  processing  problems  requiring 
enormous  amounts  of  parallel  data  input  to  the 
processor  but  cos^aratively  little  data  flow  out. 

An  excellent  example  of  this  type  of  processor  is 
the  systolic  array  described  in  the  previous 
section.  Let  the  input  and  output  vectors  be  of 
length  N,  and  the  matrix  be  NxN.  The  elements  of 
the  input  vector  can  be  entered  over  a  serial  elec¬ 
tronic  communication  line  and  likewise  the  eleaenta 
of  the  output  vector  can  leave  the  chip  via  a  single 
serial  communication  line.  However,  the  (2N-1) 
processors  mutt  receive  an  average  of  H  different 
matrix  elements  every  clock  cycle,  thus  requiring 
either  an  enormous  number  of  parallel  input  channels 
or  •  large  degree  of  multiplexing,  tdtich  will  slow 
down  the  speed  of  operation.  •  To  provide  parallel 
electronic  input  channels  requires  that  the  chip 
have  a  very  large  nuaber  of  pins.  This  rather 
severe  price  can  be  avoided  if  data  is  input  via  a 
corresponding  mnaber  of  imaging  optical  inter¬ 
connections,  as  shown  in  fig.  20.  It  is  only 
necessary  to  integrate  detectors  on  the  chip.  The 
need  for  a  large  nuaber  of  bonding  pads  is  aliain- 
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(t«4,  VlhiU  parallel  electronics  are  needed  to 
drive  the  light  source*  Chat  are  iaaged  onto  the 
chip,  such  electronics  can  be  relatively 
aacroscopic  in  site,  free  from  the  constraints 
present  in  a  single-chip  environaent.  The  priae 
computational  eleasnts  of  the  architecture  remain 
in  the  single  chip  environaent,  as  needed  for  high 
speed  and  low  power  dissipation. 


OPTICAL  INPUT  FOR 
MATRIX  ELEMENTS 
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Figure  20 


SYSTOLIC  PROCESSOR 


Replacement  of  pins  by  iaaging 
optical  connections. 


The  above  ideas  are  clearly  speculative  but  they  do 
provide  one  possible  scenario  for  future  directions 
in  optical  data  processing.  Here  it  is  no  longer 
optics  that  does  the  actual  computing  but  the  role 
of  optics  as  a  means  for  interconnections  is  no 
less  important  than  that  of  the  digital  electronic 
processors  that  perform  the  nuaerical  computations. 
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ABSTRACT 

The  statistical  properties  of  the  spatial  derivatives  of  the  phase  of  a  mono¬ 
chromatic  speckle  pattern  are  studied.  Initially,  a  one-dimensional  probability 
density  function  for  the  derivative  of  the  phase  is  obtained  and  compared  to  the 
solution  for  the  analogous  problem  concerning  instantaneous  frequency  of  nar¬ 
rowband  Gaussian  noise.  Subsequently!  *  two-dimensional  probability  density 
function  b  derived  which  depends  on  the  two  first  and  three  second  spatial 
moments  of  the  illumination  intensity  distribution  of  the  scattering  object.  Some 
sample  intensity  distributions  are  considered  for  which  explicit  expressions  for  the 
probability  density  function  are  given. 
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Introduction 

Speckle  patterns  arise  when  highly  coherent  light  b  transmitted  through,  or  reflected  from, 
an  object  with  a  surface  that  b  rough  on  the  scale  of  a  wavelength.  It  has  been  found  fruitful  to 
treat  the  fields  present  in  such  a  pattern  as  a  random  process,  the  randomness  being  over  an 
ensemble  of  macroscopic  ally  similar  bat  microscopically  different  rough  surfaces.  In  thb  paper  we 
are  concerned  with  a  very  specific  statistical  property  of  such  patterns,  namely  the  statistical  dis¬ 
tribution  of  geometrical  ray  directions  contained  in  the  scattered  fields.  Equivalently,  we  are 
interested  in  the  joint  probability  density  function  of  the  two  components  of  the  gradient  of  the 
phase  dbtribution  in  a  speckle  pattern. 

The  statbtical  dbtribution  of  ray  directions  b  of  interest  in  a  number  of  practical  problems. 
For  example,  when  a  hologram  b  formed  by  interference  of  a  plane  reference  wave  and  a  wave 
generated  by  a  diffuse  object,  knowledge  of  the  statbtical  dbtribution  of  phase  slopes  allows 
specification  of  the  statbtical  dbtribution  of  local  spatial  frequencies  on  the  recording  medium.  If 
the  hologram  b  thick,  then  thb  information  also  allows  specification  of  the  local  tQt  of  the  fringes 
within  the  emulsion,  a  quantity  that  influences  the  average  diffraction  efficiency  that  can  be 
obtained.  Finally,  the  results  may  be  of  interest  in  the  field  of  adaptive  optics,  since  many  wave- 
front  sensors  measure  the  gradient  of  the  phase  of  the  incoming  wave. 

To  our  knowledge,  the  statbtics  of  ray  directions  in  speckle  patterns  have  not  been  dis¬ 
cussed  in  the  literature  before.  Neither  of  the  two  general  surveys  of  thb  field  discusses  the  prob¬ 
lem  |l,2].  The  most  closely  related  works  are  those  of  Ebeling  (3)  ,  who  investigated  the 
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statistical  properties  of  the  spatial  derivative  (along  one  arbitrary  direction)  of  the  amplitude  and 
intensity  in  a  speckle  pattern,  and  Ohtsubo  (4}  ,  who  studied  the  zero-crossing  rate  of  the  deriva¬ 
tive  of  the  intensity  of  a  speckle  pattern.  Also  worth  mention  are  the  classic  works  of  Longuet- 
Higgins  (5]  ,  who  studied  the  statistical  properties  of  the  derivatives  of  wavefronts  with  Gaussian 
phase  deformations  but  constant  intensities.  Our  problem  is  significantly  different  in  that  both 
the  intensity  and  the  phase  are  random,  and  the  phase  statistics  are  not  Gaussian.  Finally,  we 
should  mention  that  our  problem  has  an  exact  one-dimensional  analog  that  has  been  studied  by 
communication  theorists,  namely  the  statistical  distribution  of  the  instantaneous  frequency  of  nar¬ 
rowband  Gaussian  noise  (see,  for  example,  reference  (6)).  While  the  methods  of  solution  of  the 
one-  and  two-dimensional  problems  are  similar,  the  two-dimensional  solution  is  more  complex, 
and  the  solution  in  two  dimensions  cannot  be  deduced  simply  from  knowledge  of  the  result  in  one 
dimension. 

Background 

A  monochromatic  speckle  pattern  can  be  described  as  a  sum  of  contributions  from  N 
independently-phased,  coherent  radiators.  A  single  component  of  one  polarization  component  of 
the  electric  field  E  at  time  t  has  the  form 

E(i,y,z,t)=-^£\ok  |  e‘V(—4  (1) 

where  |  a*  |  and  d*  *re  the  amplitude  and  phase  of  the  hth  radiator  and  are  functions  of  x,  y, 
and  x.  6  represents  the  initial  phase  and  b  a  random  variable  uniformly  distributed  from  -a  to  x. 

Separating  E  into  real  and  imaginary  parts,  we  have 

E(x,p,s,t)  —  r  +  i  (2) 

1  s 

■*  ”7=r  2  I fl*  I  |(cosd*cos(wf  +  0)-sind*sin(u/l  +  0))  +  t(cosd*sin(uf  +  $)+  sind*cos(u;t  +  0))\ 
vrv 

-«  J  o  |  |cos(wf  +  d)  +  isin(wf  +  d)J 

where  d  is  the  phase  of  the  resultant  wave.  Statbtical  properties  of  speckle  are  studied  using  the 
following  assumptions: 
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Consider  r[xi,y,x,t),  i(z„y,z>l),  r[x2,y,x,t),  and  i{zj,y,z,t).  The  central  limit  theorem 
tells  us  they  are  jointly  Gaussian  random  variables  (7).  It  is  also  true  that  r(z,y,z,t), 
3  3 

—  r(x,y,:,t),  i(x,y,x,t),  and  —  i(z,y,z,f)  are  jointly  Gaussian  [8,  p.475j.  As  a  shorthand 

notation,  ~^f(*,y,z,t )  will  be  written  as  f,(z,y,z,t).  The  general  form  of  the  distribution  den¬ 
sity  function  is  [8,  p.  255] 


where 


and 


p(r,*,r„it) 


1 

4*V3et H 


l“l  “  I  r  i 


r.  41 


Af 


<rr>  <ri>  <rrM>  <ri,> 

<»>>  <»»>  <»Y,>  <«,> 

<rtr>  <r,f>  <r„r,>  <r,i,> 
<*,r>  <i,i>  <i,r,>  <«,*•> 
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Following  Ebeling  [3]  ,  we  note  two  relations  in  order  to  evaluate  elements  of  M : 


(0 


<A{zx,y,x)A{x2,y,z)><tiU>tt%*  (10) 

=  0 

since  <eiU>  «  0.  Note  A(x,y,z)  is  the  complex  amplitude  of  the  electric  field  given  in 
Eq.  (1).  Therefore,  equating  real  and  imaginary  parts,  we  find  that 

R„{*x)  -  tf„(A*)  (II) 

Ari(Ax)  —  -J?<,(Az)  • 

nfi 

(2)  From  Papoulis  (8,  p.317| ,  with  representing  yj-/  ,  we  have 

“  (“O'"  m)  •  02) 
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M  is  symmetric  and  the  upper  2x2  portion  of  it  b  already  known,  so  seven  elements  remain  to  be 
determined.  Note  that 


Similarly, 

| 

<*Tsrt^  "  -  1  A»~0  B  K  • 

(13) 

<*.*«>  -  *.  • 

From  Eq.  (11),  it  b  seen  that 

(M) 

<»»*.>  ' 

and  thus 

^  1  t?  * 

™  —  ^(^*)  1  “  ^p,Afr(Ai)  J  4x«o  “ 

(15) 

Continuing,  since 

<r,i,>  ==  <»,ra>  wp  . 

we  find 

<//>»  _  <(//)*>  =  2<//»> 

<rr->  "  i  to<"> 

”  i  ^(co“u“t)“0  ■ 

(16) 

By  the  same  reasoning,  <  ii,  >  —  0  as  well  Finally, 

<n,>  =  -^^„(Ax)  |  =  c. 

(17) 

and  <r,i>  b  simply 

-<rig>.  Therefore, 

0s  0  0  e, 

0  o*  -e,  0 

M  -  (18) 

0  -c,  la  0 
c.  0  0 

»  t 

From  thb  form,  the  determinant  and  inverse  of  M  are  found  to  be 

dei  A/  -  (o^.-c.Y  (19) 


•  6- 


and 


Ml  -  ~ 


4,  0  0  -c, 

0  4,  e,  0 

0  e,  «*  0 

|— c,  0  0  o* 


(20) 


where 


D  =  V3et  A/  =  o*4,  -  e*  . 

Performing  the  matrix-vector  multiplications  of  Eq.  (7),  we  get 

P{',i,rt,h)  =  ^f^exp{-^(4,(r*+  »*)  +  o*(r  *+  »*)  +  2c,(»r,-r^)J} 

The  transformation  of  variables  to  intensity  /  and  phase  ^  b 

r  —  \/7  cos  (wf +  d) 

i  we  VI  *in(wl+  4) 

r,  =  -^^cos(w<+  4)  -  VI  *,sin(wf+  $) 

*,  “  ~^pm[ut+  4)  +  VI 4,cos(ut+  d) 

IUII  -  -j 

with  a  resulting  probability  density  function 


(21) 


(22) 


(23) 


*(/.*,/.,*.)  -  «P{-^jl4.J  +  o*(~+/dJ)-2c,/d,l}  . 


(24) 


Since  the  probability  density  function  of  d,  b  our  goal,  H  b  necessary  to  integrate  over  /,  d  and 
/,.  Thb  b  a  straightforward  process  and  results  in 

(g»4,  -  e*) 


>(*.)  -  7T 


where,  to  repeat, 


2a  (o*tf-2c,+,+  |,JV* 


e*  —  0) 

*•  "  “  ftrr(A*)  I  Ai-S 


(25) 


(26) 


-7- 


This  expression  cm  be  compared  with  the  distribution  obtained  in  the  analogous  problem  con¬ 
cerning  the  instantaneous  frequency  of  narrowband  Gaussian  noise.  Blachman’s  formula  [8,  p.6l] 
is 


f(*») 


2 xV_ 


(*.*  +  4  Wf* 
where  p*  b  the  spectral  width  and  is  defined  as 


,  1  <r,ra> 

Therefore,  the  relationship  to  our  notation  b 


(27) 

(28) 


*.  —  4x/V  (28) 

In  addition,  in  Blachman’s  formulation  of  the  problem,  c,»0  because  he  tales  as  hb  reference 

frequency  the  centroid  of  the  noise  spectrum.  Hence,  it  b  seen  that  the  two  formulas  are 

equivalent. 


Probability  Density  Function  of  Slopes  In  Two  Dimensions 

Our  ultimate  goal  b  to  solve  the  two-dimensional  problem,  finding  an  expression  for 
p(4i>4f)-  Thb  will  give  us  the  statistical  distribution  of  geometrical  ray  directions  since,  if  a  and 
p  are  the  direction  cosines  with  respect  to  the  *-  and  y-axes,  then 


Extrapolating  from  the  previous  section,  the  desired  density  function  b  seen  to  be  obtainable  from 
y(r»*«r»«'*<r#»V)'  earlier  reasoning,  these  six  random  variables  are  jointly  Gaussian.  Hence, 
the  general  form  of  the  dbtribution  b 

p(r,f,r„f„r„i,)-  exp(-j(u|Af-,(u|r)  (31) 

- 


(32) 


and 


M 


\ 

lul  “  l r 

‘  f, 

r» 

*»  1 

<rr> 

<ri> 

<rr,> 

<r»j> 

<rrt> 

<ri,> 

<»>> 

<ii> 

<«>,> 

<«»,> 

<ir,> 

<»»,> 

<r,r> 

<r,i> 

<r,r,> 

<»■.».> 

<r,rt> 

<r.S> 

<»,r> 

<*,»> 

<*,*,> 

<*A> 

<r,r> 

<r,i> 

<r,r,> 

A 

V 

<r»S> 

<ifT> 

<•;»> 

<•>.> 

A 

V 

<*»r»> 

<vs> 

(33) 


Because  of  symmetry,  twenty-one  elements  most  be  found,  of  which  ten  have  already  been  deter¬ 
mined.  The  remainder  of  the  elements  are  discussed  in  Appendix  A.  The  resulting  covariance 
matrix  has  the  form 

jo*  0  0  c,  0  e, 

0  o*  -t,  0  —Cg  0 

0  -e,  4,  0  dO 

e,  0  0  I,  0  d 

0  -c,  d  0  ft,  0 

c,  0  0  d  0  *, 


M 


(34) 


where 


d  ■  <r,rf  >  mm  - 


0  d 


flA*  d*tR,,{AlAy'1 1  &f 


(35) 


The  inverse  of  Af  is 


I 


D  -  V$*H  **0*b,i,-o*d*-i,c*-itc*+2c,c,4  ..  (37) 

Substituting  into  Eq.  (31),  we  obtain 

P(r,»\r„»,,r„if)  -  »*)  +  (a2i,-c/Xr,*+  if)  +  (38) 

{afi,-cf)(r*+  if)  +  2 {c,d-k,c,)(rit-ir,)  +  2(c,i-l,  )  +  2(c,«*-o*dXf«'»+ ».»,))} 

Again,  a  tranafonnation  is  seeded  to  get  the  distribution  of  /,  4>,  /„  ,  and  .  Tbe  rela¬ 

tions  are : 

r-/T  cos(w<+ 

»  ■  VTria(ut+  4) 

r ,  “■  cos(wl+  4) -VI d,sin(wt+  $) 

•«  —  -^jtia(ut+  4)+  VI  d,cos(wl+  <i)  (39) 

r,  —  ^  ’<«(wl+  *)-  \/7 *ftin(wH-  <>) 

^  —  •^=-sin(wJ+  4)  +  VI ^,eos(w<+  4)  . 

The  following  transform  atioa  expressions  result : 

r*+i*-/ 


10 


ri,  -  ir, 
ri,  -  ir, 

r,rt  +  »,V 

«/ii 


i*. 

Ml. 

4 1 
_1_ 

8 


+  I***, 


(40) 


The  probability  density  function  hence  becomes 

(o*‘,-«,*K^+ /#.*)+  (41) 

Integration  over  /,  *,  /»  and  /,  is  performed  in  Appendix  B.  The  resulting  joint  probability  den* 
sity  function  of  the  e-derivative  of  the  phase  and  the  y -derivative  of  the  phase  b 

Pi*..*,)  -  ~-\{b,b,-d?)  +  (42) 

+  2(*<M,c,)<>,  +  2 M-*,e,)*,  +  2(c,c,Vd)*,d,)-8  . 


Examples 

It  b  instructive  to  calculate  p  (*,,*,)  for  some  basic  htensity  distributions  at  the  scattering 
plane.  The  parameters  of  the  dbtribution  in  Eq.  (42)  depend  on  spatial  correlation  functions.  It 
b  possible  to  rewrite  them  all  as  depending  on  the  autocorrelation  function  of  the  complex  ampli¬ 
tude  /^.(Az.Ay).  Thb  can  be  seen  by  first  noting  that 

Ra.(&xAv)  -  /?*,.( Az.Ay)  (43) 

By  expanding  E  in  terms  of  its  real  and  imaginary  parts,  the  following  relations  result : 

J?„(Az,Ay)  —  yRe  tf^.fAz.Ay)  (44) 

/frt(Az,Ay)  —  -|lm  /^.(Az.Ay)  . 


«*—  ■jRe/fiM.(Az— 0,Ay— 0) 

*.  -  -  [flAz*  1 *;-•] 


Hence, 
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'■  “  71"  I  £-] 

*'  -  -iR‘  1  s;=s]  <“> 

I  £-1 

1  “  ~lKt  I  ijs]  ‘ 

The  last  relation  needed  to  calculate  p  {$,,$,)  fa  provided  by  the  Van  Cittert-Zernike 
theorem  (0|, 

R^sAv)  -  rhFT,{l{u,v))M  (40) 

**  Xi’U 

where  FT 2  denotes  2-D  Fourier  transformation,  /(«,«)  fa  the  illumination  intensity  distribution  at 
the  scattering  plane,  and  r  is  the  distance  between  the  scattering  plane  and  the  observation  plane. 
Combining  the  results  of  equations  (45)  and  (40),  we  And 

*■  “ 

—  -rr//«/(u,»)d«dv  (47) 

as  e 

cr  -  •jy //«*(«.•)*•* 

d  -  ^n-/7«v/(.,v)**  * 

AS  • 

As  a  first  example,  consider  a  rectangular  spot  of  dimension  LxW  .  The  intensity  distribu¬ 
tion  is  written  as 

/(u,v)-  a  r«ct|jjrect|~j  (48) 

where  «  is  a  constant  and  rect  (/)  is  defined  to  be 
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ttct(f) 


1  for|/j<i 


0  *°r|/|>7 


Taking  the  Fourier  transform,  we  obtain 


R^xAv)  “  «  ^7  ,inc(iAx)si,,c(^A») 


where 


sine  (/)< 


ESllIl 


It  is  convenient  to  define  a  “speckle  width",  in  each  direction,  as  the  distance  in  which  the  corre¬ 
lation  function  falls  to  tero.  Let 


6x  ■=■  speckle  width  in  x -direction  =■  — 

L 

6y  =  speckle  width  in  y-direction  =»  —  . 


The  parameters  needed  for  Eq.  (42)  are  determined  from  this  equation  to  be 


Hence, 


2 v*m 

»  « — ««?— . 
’  W*ftb) 

c.-0 

*  . — 

e,  -0 

d  -0  . 


"  3{*r)W  [(dp)*  + 


W’fthf, 
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If  L*"W  (and  therefore  6*  ■  fy) , 


3  [6*?\ 


A  graph  of  this  function  is  shown  in  Fig.  1.  for  the  following  parameters: 

X  —  .5  •  l<r*  m 
i>ln 

L  —  1  am  (50) 

Sx  mm  .0005  m 

Both  snd  mage  from  -10,000  to  10,000  radians  per  meter  with  the  (0,0)  point  at  the  center 
of  the  graph. 

As  a  second  example,  a  symmetric  Gaussian  scattering  spot  is  considered.  The  intensity  dis¬ 
tribution. 


/(«.*)  —  «  e  l'(*  , 

b  defined  so  that  Lj 2  is  the  l/e  point  in  each  direction.  The  autocorrelation  function  is 
R^CxtAv)  -  «*  exp|-*»(Ax*+  • 

Let  the  definition  of  speckle  width  in  this  case  be  the  l/e  point.  Thus 


and  the  nonsero  parameters  of  p(d*>d, )  ■**  found  to  be 

ir(bf 

The  resulting  distribution  density  is  * 


14. 


Conclusion 

The  joist  probability  density  function  of  the  *-  and  y-  components  of  the  phase  in  a  mono* 
chromatic,  fully  developed  speckle  pattern  has  been  found.  Equivalently,  the  statistical  diatribe* 
tion  of  ray  direction?  in  suck  a  speckle  pattern  have  been  derived.  The  results  show  that  the 
statistics  in  question  depend  os  only  the  two  first  and  three  second  spatial  moments  of  the  illumi¬ 
nation  distribution  of  the  scattering  object  but  not  on  other  details  of  that  distribution.  This 
problem  is  the  two-dimensional  analog  of  the  problem  of  finding  the  statistics  of  the  instantaneous 
frequency  of  narrowband  Gaussian  random  noise.  The  results  should  prove  useful  in  studies  of 
the  average  diffraction  efficiency  of  thick  holograms  of  diffuse  objects  and  may  be  of  interest  in 
adaptive  optics. 


Acknowledgments 

The  authors  gratefully  acknowledge  support  from  AFOSR  and  an  IBM  predoc tora!  fellow¬ 
ship. 


Appendix  A 

The  6x6  covariance  matrix  M  (Eq.  (33)  in  the  main  text)  has,  on  account  of  symmetry, 
twenty-one  distinct  elements  which  must  be  found.  Ten  of  these  were  discussed  In  the  text. 
Furthermore,  seven  more  have  exact  analogies  to  cases  previously  studied.  Thus  : 

<rs r»>  -  -7£"s*"(A*)  1 0,-0  «  (61) 

and  <*f «,  >  »  I,  m  well; 


<r,v>-0 
<rrt>  -  <»i,>  —  0 

<*»>  “  I  Ar-s  ■  ** 

and  <r,i>  —  -<ri#>. 


(62) 

(63) 

(64) 
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To  study  the  remaining  four  covariances,  it  b  necessary  to  find  a  relation  similar  to  Eq. 
(12),  but  where  the  partial  derivatives  of  the  two  functions  to  be  cross-correlated  are  with  respect 
to  two  different  variables. 


Following  the  reasoning  of  Papoulis  [8,  p.316],  we  can  evaluate 

—  E{!t[xl,yl)gt{z3,y2)) 

in  two  separate  steps.  First,  we  observe  that 


(85) 


f  (»«,»!+  thtjzj'Vi)  j|  */,(*  i,*7.Vi.Vi+  tH?/«(*r.*».yi.lb) 


Taking  the  limit  of  both  sides  as  <-*0,  we  get 


o 


Second, 


(66) 


(67) 


A 

4 

i 


f'{,^  gar  .  (6g) 


Again  taking  the  limit  as  e-»0,  we  find 


d  0 


6*i  8yt 

If  the  processes  are  jointly  spatially  stationary,  then 


where 


*/,#,<*».**)  - -  99-x  ^i?/,(A*. Ay) 


A*  —  *»  -  ** 
Af  —  Vi-ti 


6 Ax  8 Ay 


8 Ax  OAm 


(69) 


(70) 


(71) 


Returning  to  the  remaining  covariance  terms,  we  have 


>16 


and  so 


Lastly, 


<r,*,>=s<  V»>  53  0  • 


<r,rt>  =  -■ 


d  d 


•J?,,(A*,Ay)  | 


AtmO  ™*  i 
Ai«e 


(74) 


flAx  5Ay 

and  <»,»',  >  «=»  <r,r(  >  .  The  resulting  covariance  matrix  b  given  in  Eq.  (34)  of  the  main  text. 


Appendix  B 

The  probability  density  function  obtained  when  the  transformation  to  six  variables  related 
to  intensity  and  phase  b  completed  b,  as  given  in  Eq.  (41)  of  the  main  text, 

P (l, ~  +  (oVc/X-£+  W)  + 

M.-c ;i£+I*;)+  2{c,d~b,c,)I$,+  2(ctd-b,c,)It,  +  2(c,c,-^rf)(^+ /*,*,)]>  . 

Integration  over  ^  yields  a  factor  of  2x .  The  integration  over  I,  b  of  the  form 

i! 


-CO  V  • 


where 


*,  q3t»~<v  »  (c,ct-o*rf)Jf 

*  "  8/2)  AID 

The  form  of  the  integration  over  /f  b  then  simply 

i  V  •  »l.  « 


It  b  useful  to  note  that 


so  a  can  be  rewritten  as 


(o*»,-c*)(o*6#-c/)-  (c,c,-o*df  —  0*1) 


e* 


(75) 


(76) 


(77) 


8/(o*»,-c/)  ‘ 


(78) 


These  integrations  leave 


p(l  =■  +  (e**.-e*)d*  +  P«) 

2(etd-htet)<t>t  +  2(e,d-J,cf  )^f  +  2(e,c,-e*d)dsd,|)  . 

The  last  integration  takes  the  form 

fle-'df  -  .  (80) 

Finally,  the  joint  probability  density  function  of  the  x -derivative  of  the  phase  and  the  >- 
derivative  of  the  phase  is 

r{K*t)  —  ~-\ (M,-*2)  +  («***-</)#  +  (e**,-c»*)tf  + 

2(ct  d-br  c,)fiM  +  2 (e,rf-i,cf)£#  +  2(c»c^-o*d)d»dfl_*  • 
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Figure  Captions 


Fig.  i  :  Joint  probability  density  function  of  the  x-  and  y-  derivatives 
of  the  phase  for  a  uniform,  square  intensity  spot  in  the 
scattering  plane. 


